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ABSTRACT 



Turbulent magnetic fields are ubiquitous in space physics and astrophysics. The 
influence of magnetic turbulence on the motions of charged particles contains the 
essential physics of the transport and acceleration of energetic charged particles in 
the heliosphere, which is to be explored in this thesis. After a brief introduction 
on the energetic charged particles and magnetic fields in the heliosphere, the rest 
of this dissertation focuses on three specific topics: 1. the transport of energetic 
charged particles in the inner heliosphere, 2. the acceleration of ions at collisionless 
shocks, and 3. the acceleration of electrons at collisionless shocks. We utilize various 
numerical techniques to study these topics. In Chapter 2 we study the propagation 
of charged particles in turbulent magnetic fields similar to the propagation of solar 
energetic particles in the inner heliosphere. The trajectories of energetic charged 
particles in the turbulent magnetic field are numerically integrated. The turbulence 
model includes a Kolmogorov-like magnetic field power spectrum containing a broad 
range of scales from those that lead to large-scale field-line random walk to small 
scales leading to resonant pitch-angle scattering of energetic particles. We show that 
small-scale variations in particle intensities (the so-called "dropouts") and velocity 
dispersions observed by spacecraft can be reproduced using this method. Our study 
gives a new constraint on the error of "onset analysis", which is a technique com- 
monly used to infer information about the initial release of energetic particles. We 
also find that the dropouts are rarely produced in the simulations using the so-called 
"two-component" magnetic turbulence model (IMatthaeus et all Il990l ). The result 
questions the validity of this model in studying particle transport. In the first part 
of Chapter 3 we study the acceleration of ions in the existence of turbulent mag- 
netic fields. We use 3-D self-consistent hybrid simulations (kinetic ions and fluid 
electrons) to investigate the acceleration of low-energy particles (often termed as 
"injection problem") at parallel shocks. We find that the accelerated particles al- 
ways gain the first amount of energy by reflection and acceleration at the shock layer. 
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The protons can move off their original field lines in the 3-D electric and magnetic 
fields. The results are consistent with the acceleration mechanism found in previous 
1-D and 2-D simulations. In the second part of Chapter 3, we use a stochastic inte- 
gration method to study diffusive shock acceleration in the existence of large-scale 
magnetic variations. We show that the 1-D steady state solution of diffusive shock 
acceleration can be significantly modified in this situation. The results suggest that 
the observations of anomalous cosmic rays by Voyager spacecraft can be explained 
by a 2-D shock that includes the large-scale magnetic field variations. In Chapter 
4 we study electron acceleration at a shock passing into a turbulent magnetic field 
by using a combination of hybrid simulations and test-particle electron simulations. 
We find that the acceleration of electrons is greatly enhanced by including the ef- 
fect of large-scale magnetic turbulence. Since the electrons mainly follow along the 
magnetic lines of force, the large-scale braiding of field lines in space allows the 
fast-moving electrons interacting with the shock front multiple times. Ripples in 
the shock front occurring at various scales also contribute to the acceleration by 
mirroring the electrons. Our calculation shows that this process favors electron ac- 
celeration at perpendicular shocks. We discuss the application of this process in 
interplanetary shocks and flare termination shocks. We also discuss the implication 
of this study to solar energetic particles (SEPs) by comparing the acceleration of 
electrons with that of protons. The intensity correlation of electrons and ions in SEP 
events indicates that perpendicular or quasi-perpendicular shocks play an important 
role in accelerating charged particles. In Chapter 5 we summarize the results of this 
thesis and discuss possible future work. 
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CHAPTER 1 



Introduction and Background 



1.1 Overview 

The heliosphere (Figure 11.11) is struc t ured by plasma flows and magnetic fields. As 



the supersonic "solar wind" (jParkerl . Il958l ) expands from the solar atmosphere, it 



forces plasmas and magnetic fields outward, forming the cavity that interacts with 
the surrounding interstellar medium. At around 100 AU, the solar wind suddenly 
slows down and forms the termination shock. The heliosphere is a natural laboratory 
for physical processes involving charged particles and changing magnetic field. The 
related physical quantities are continuously measured by spacecraft and ground- 
based instruments. In particular, the Voyager spacecraft are currently exploring the 
heliopause, which is the boundary between the heliosphere and the local interstellar 
medium. It marks the frontier of the solar system and the farthest distance (~ 121 
AU at the time of this writing) that man-made instruments have ever reached. 

Energetic charged particles are a minor component of space plasmas, but have 
important and profound effects. They are serious concerns to space weather because 
of their hazardous effects to astronauts and space satellites. Energetic charged 
particles propagate at high speeds and carry important information about their 
source regions and the media they propagate through. Since the early measurements 
by Victor Hess in 1912, scientists have been measuring energetic particles with 
energies up to about 10 21 eV and the electromagnetic radiation they produce. 

The acceleration and transport of energetic charged particles are fundamental 
problems in space physics and astrophysics, in which electric and magnetic fields 
often play an essential role. From the first principle, the motions of charged par- 
ticles in electromagnetic fields are governed by the Lorentz force. The large-scale 
transport and acceleration of energetic charged particles are often described by the 
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Figure 1.1: The interaction between interstellar medium and the heliosphere. Figure 
provided courtesy of Steve T. Suess. 
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Parker's transport equation (jParkerl . Il965l . see Equation 11.131) . and the effect of 
turbulent magnetic fields is approximated by a large-scale spatial diffusion tensor 



( I Jokipii . 



1966 



197ll ). However, recent studies suggest that the simple spatial dif- 



fusion approximation may be oversimplified and the observed behavior of energetic 
particles is often more complicated. For example, during some solar energetic par- 
ticle events, the intensity-time profiles often show small-scale sharp variations. In 
addition, when the velocity distribution of charged particles is highly anisotropic, 
it cannot be properly described by the transport equation that assumes a quasi- 
isotropic distribution. It is important to adequately describe this for problems such 
as the evolution of the velocity distribution of energetic particles during their propa- 
gation and the acceleration of low-energy particles at shocks. Understanding this is 
a challenge to theoretical studies due to the complex nature of particle trajectories 
in turbulent magnetic fields. Thanks to recent advances in computing capabilities, 
the acceleration and transport of charged particles can be studied by numerical sim- 
ulations. This dissertation focuses on understanding the transport and acceleration 
of energetic charged particles in the existence of turbulent magnetic fields. 



1.2 Charged Particles in the Heliosphere 

The heliosphere is permeated by charged particles of different origins. There are 
two important components of charged particles that contribute to the global dy- 
namics inside the heliosphere, i.e., the solar wind and pickup ion s. The solar w ind 
is a continuous plasma flow coming from the solar atmosphere (jParkerl . Il958l ). It 
is accelerated to supersonic speeds close to the Sun and propagates outward. The 
solar wind is often divided into two distinct components, termed as the slow so- 
lar wind and fast solar wind. The fast solar wind represents a plasma flow with 
a temperature about 8 x 10 5 K and a speed of about 800 km/s, whereas the slow 
solar wind represents a hotter ( T ~ 1.5 x 10 6 K) a nd denser plasma flow with a 
slower speed of about 400 km/s ( Meyer- Vernetl . 120071 ). The observations of the solar 
wind at different latitudes by the Ulysses spacecraft have shown that the slow solar 
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wind is confined to "streamer belt" that is about 20 degrees around the heliospheric 
current sheet at solar minimum. At the same time the fast solar wind entirely oc- 
cupies higher lati tudes. At solar maxim um, the solar wind becomes more mixed 



and complicated ( iMcComas et al. 



inter stellar neutral particles (|Axford 



2003|) 



1972 



ickup ions are mainly origin ated from 



Vasvliunas fc Siscod. 



1976 



Gloeckler 



et al., 



1995 



19951). with mino r contributions from inner sou rce pickup ions ( IGeiss et al. 



Gloeckler &: Geissl. 



( Ilpavich et al. 



1986 



1998 



Schwadron et al. 



Gloeckler et al. 



20001 ) and pickup ions from comets 



19861 ). Interstellar neutral particles can freely 



penetrate into the heliosphere before they are ionized by charge exchange and/or 
solar ultraviolet radiation. Once the neutral particles are ionized, they are influ- 
enced by electric and magnetic fields (see Section [L3]) . As illustrated in Figure [TT2| 
when the magnetic field vector has a component that is perpendicular to the solar 
wind velocity vector, the electric and magnetic fields embedded in the solar wind 
force them to accelerate and make gyro-motions in the frame co-moving with the 
solar wind. This is the so called "pickup" process. In observer's frame, the pickup 
ions have velocities from zero to two times of solar wind speed. After the pickup 
process, the gyroradii of pickup ions are several times larger than that of the solar 
wind particles. The gyro-motion forms a ring velocity distribution perpendicular 
to the ambient magnetic field. The distribution is unstable and expected to gen- 
erate electromagnetic waves. The "pickup ions" will be sca t tered into an isotropic 



distribution by the electromagnetic waves ( iWu fc Davidson 



1972 



Wu et al. 



1986|) 



and/or background turbulence. Starting from a heliocentric distance of about 30 
AU, the pickup ions play a dominant role in the physics of the outer hel iosphere by 
dominating the pressure of the plasma flow (IRichardson fc Stone! . 120091 ). It should 



be noted that neutral particles could a 



which has been discussed extensively (IZank 



so influence the dynamics in the heliosphere, 



1999 



and references therein). 



Energetic charged particles are a high-energy, non-thermal component of the 
plasmas in the heliosphere. They carry large kinetic energies, so their speeds are 
much faster than background fluid. They have significant effects on space weather 
and their physics is important to consider. When the energy density of energetic 
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charged pa rticles is large enough, they may even mediate the dynamics of plasma 
flows (e.g., iFlorinski et al.l . 120091 ). Energetic charged particles in the heliosphere 



have different components depending on their acceleration sites, energy ranges, 
charge states, and elemental compositions, etc. Figure [L3l illustrates various types 
of charged particles in the heliosphere, their acceleration regions, and their typical 
energy spectra. In this figure, the solar wind with energies in keV range repre- 
sents the kinetic energy of the background flow in the heliosphere. Solar energetic 
particles (SEPs) are usually observed to be from several hundred keV/nucleon to 
tens of MeV/nucleon for typical events, and even more than 1 GeV/nucleon in ex- 
treme events. The SEP events are usually divided into two classes based on the 
progress made i n the las t thre e decades, termed as "impulsive events" and "grad- 

and references therein). The SEPs associated with 



ual events" (see 



Reames 



1999 



impulsive events are thought to be accelerated in solar flares. Impulsive events are 
characterized by the impulsive peaks in their intensity-time profiles, confined source 
regions in longitude, high ionization states, and overabundance in isotope ratios 
such as 3 He/ 4 He and Fe/O compared with the values in the solar wind. Energetic 
particles related to gradual events are thought to be accelerated by collisionless 
shocks driven by coronal mass ejections (CMEs). The gradual events have extended 
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intensity-time profiles because of the continuous acceleration at shock fronts. They 
are also more w idely distributed in longitude, indicating a broadened source region 
( IReamesl . Il999l ). Recently, this bi- model picture has been challenged. It has been 
found that for most events, SEPs appear to have a mixed property of the two classes 
of events. For example , some large gradu al SEP events show a substantially high 
ionization c harge state (Mazur et al.L Il999h and enha nced isotope ratios in 3 He/ 4 He 



and Fe/O (I Cohen et al. 



1999 



Mason et al. 



19991 ). In large SEP events, flares 



and CMEs usually occur together, therefore one may expect that both of the pro- 
cesses play a role. Identifying their relations and contributions to large SEP events 
is still under hot debate. Energetic storm particles (ESPs) are usually associated 
with the passage of travelling interplanetary shocks, where the ions can often be 
accelerated to from several tens of keV/nucleon to MeV/nucleon, and occasionally 
more than 10 MeV/nucleon. Energetic charged particles associated with corotating 
interaction regions (CIRs) formed by the interaction between the fast and slow solar 
wind streams are often observed. These particles appear to have a higher energy 
range compared to that of ESP s. Anomalous co smic rays (ACRs) are thought to 
be originated from pickup ions (IFisk et all Il974l ). and are accelerated to energies 
betw een seyeral MeV/nucleon to 100 MeV/nucleon at the termination shock f Pesses 



et al., Il98ll ). Galactic cosmic rays (GCRs) that have energies typically larger than 
10 MeV/nucleon are coming from the outside of the heliosphere. 

One remarkable feature in Figure 11.31 is that the energy spectra of the acceler- 
ated ions are often close to a power law, which indicates a common acceleration 
process. The acceleration of electrons is also observed, and sometimes accompanies 
the acceleration of ions. We will discuss the acceleration of electrons in Chapter 4. 
It is worthwhile to note that many, if not all, of the energetic charged particles are 
associated with shock waves. For example, it is now established that m any SEP 
event s, especially gradual events are associated with CME-driven shocks (IReamesl . 
19991 ). Energetic particles in impulsive events are accelerated in solar flares. The 
mechanism for particle acceleration in solar flares is not clear so far. In Chapter 
4 we show that collisionless shocks are a possible candidate for the energization of 
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Energetic Particles in the Heliosphere 
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Figure 1.3: Various types of energetic charged particles in the heliosphere. The up- 
per panel indicates related physical processes in different regions of the heliosphere. 
The bottom panel indicates typical energy spectra of energetic charged particles 
obse rved at 1 AU. See text for more details. The pictures are adapted from Stone 



et al. (119981 ) with permission by Springer Science + Business Media. 
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charged particles in solar flare regions. Energetic particles can be accelerated at 
interplanetary shocks driven by interplanetary CMEs and CIRs. AC Rs are thought 



to be pickup ions accelerated at the solar wind termination shock fjPesses et al 



1 9 8 ll ) . GCRs are usually thought to be accelerated by astrophysical shocks such 
as supernova blast waves. Their spectrum is observed as the remarkable power law 
spectrum of cosmic rays. 

After the acceleration, the energetic particles travel in and interact with the he- 
liospheric magnetic field. Understanding the propagation of energetic particles is 
difficult since, in general, the motion of charged particles in a turbulent magnetic 



field is very complicated. The transport of energetic parti cles in 



netic field is usually considered to be a diffusion process (IJokipii 



l eliospheric mag- 



1966 



1971( 1. In 



principle, if the motion of energetic particles is well understood, they can be used 
as a tracer of magnetic field structure in the heliosphere. 



1.3 The Heliospheric Magnetic Field and its Fluctuations 

Since plasma flows in the heliosphere are highly electrically conductive, the helio- 
spheric magnetic field is frozen in the background fluid to a high degree. It can be 
inferred from the generalized Ohm's law that the only macroscopic electric field in 
this situation is due to the motional electric field E = -U f u„ x B/c fsee. e.g.. Krall 



& Trivelpiece, Il973l ). where Ufi ow is the background flow speed, B is magnetic field 
vector, and c is the speed of light in vacuum. The supersonic solar wind carries the 
magnetic field to many AU. Because of solar rotation, the magneti c field lines of 



1958 



Hundhausen . 



force have a spiral shape termed as the Parker spiral (iParker 
19721 ). In spherical coordinates (radial heliocentric distance f, polar angle 9, and 
azimuthal angle <j)), the average magnetic field is given by 



p2 

B = B r r + Brf = B.-l-t 



rtt s sm9 _ ig. 



(1-1) 



v ' V S w r< 
where B r and B$ are the magnetic fields in the f and direction respectively, B s 
is the radial magnetic field close to the Sun at heliocentric distance R s , Vsw is the 
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speed of the solar wind, and tt s is the angular frequency of the solar rotation. In 
this steady-state model, B s , Vsw, and Q s are assumed to be constants. 

At a low heliocentric latitude and a distance of 1 AU, the average angle be- 
tween the direction of magnetic fields and the orientation of the solar wind flow at 
low latitude is about 45 degrees. In the outer heliosphere, the plasma almost flows 
transverse to the mean magnetic field. At high latitude regions (sin^ ~ 0), the 
azimuthal component of the Parker spiral magnetic field is small. However, it is ex- 
pected that at a large heliocentric distance, the transverse perturbation of magnetic 
field lines of force dominates since it decays as B tp oc 1/r. The distant magnetic 
field almost always transverses to the radial direction while t he av erage magnetic 



field is still the Parker spiral magnetic field fjJokipii fc Kota 



19891 ). The inferred 



large scale fluctuations ha ve been observed by the Ulysses spacecraft ( IJokipii et all . 



1995 



Balogh et all Il995l ). 



Magnetic fields in the heliosphere are turbulent (ITu &; Marschl. 



1995 



Goldstein 



et al., Il995l ). It is often convenient to write the magnetic field as a summation 
of an average component and a turbulent component B = Bq + <5B. In the solar 



wind , the magnetic fluctuations are observed to be highly Alfvenic (IBelcher Sz Davisl . 



1971 



, Figure [T~4"j) . i.e., the magnetic fluctuation vector is highly correlated with the 
velocity fluctuation vector. Observations show that the power of turbulent magnetic 
field Pr versus spatial wave number k is close to a Kolmogorov law Pb oc /c~ 5 / 3 
(jColemanl . Il968l . Figure fL4|) . The power spectrum suggests that most power of the 
fluctuations is in large spatial scales, and cascades into smaller and smaller spatial 
scales until dissipation effects are important. The correlation length is observed to 
be on the or der of 10 6 km at 1 AU and increases in the outer heliosphere f Matthaeus 



& Goldstein, 



1982) • Recent theories, numerical simulations, and observations have 



revealed that the solar wind magnetic turbulence is anisotropic. In other words, 
most of the fluctuations have wave vectors transverse to the mean magnetic fi eld 



[e.g., 



Matthaeus et al. 



1990 



Goldreich fc Sridhar 



1995 



Matthaeus et al. 



1996|). 



Electric and magnetic fluctuations can also be produced by plasma instabilities. 
For example, freshly created pickup ions can have a ring distribution and excite 
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Figure 1.4: The spacecraft measurement of the Solar wind turbulence. Left: An 
exam ple of magnetic field a nd velocity measurements in the solar wind (adapted 
from iBelcher Sz Davisl . Il97ll ). Right: An example of power spectrum of magnetic 
field component measured in the solar wind fluctuations (Figure 6 in Coleman! 
Il968h . The figure is reproduced by permission of the A AS. 



ion-cyclotron waves (IWu fc Davidson 



1972 



Wu et al. 



19861 ). In the upstream 



medium of a collisionless sh ock, the streami ng of energetic particles may excite 



electromagnetic fluctuations (ILee 



1982 



19831 ). The fluctuation may be important 



for particle acceleration at quasi-parallel shocks. 



1.4 Collisionless Shocks 

Collisionless shocks are believed to be important accelerators for charged particles 
in the heliosphere. In this section we introduce the basic concepts of collisionless 
shocks. The acceleration of charged particles will be discussed in Chapter 3 and 
Chapter 4. 

Shocks are characterized by sharp transitions in the physical quantities of 
medium such as flow speed, density, magnetic field, and temperature, etc. Since 
Coulomb collisions are too infrequent in space plasmas, the kinetic energy of plasma 
flow is dissipated at shocks through other mechanisms such as the interaction be- 
tween particles and plasma waves. The shocks are termed as collisionless shocks. 
Dependent on the angle between upstream magnetic field vector and shock normal 
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Qbu-, collisionless shocks can be divided into quasi-perpendicular shocks {6 Bn > 45°) 
and quasi-parallel shocks {Obu < 45°). In the limit of ideal magnetohydrodynamics 
(MHD), shocks are a kind of discontinuities. One can relate upstream and down- 



stream medium using MHD conservation laws and Maxwell equati ons. The resu 



gives the well-known jump conditions for MHD discontinuities (e.g.. iBurgesd . Il995l ): 



t 



Bin = B 2n (1-2) 

Ui n Bi t = U 2n B2t (1-3) 

piUm = p 2 U 2n (1.4) 

pJJ^ + p. + S. = p2 ul + p 2 + ^i (1.5) 

= p 2 U 2n U 2t - ^(1.6) 



In' 





BlnBu 


47T 


(Ui 


Bl) 47T 


(U 2 


4:71 



4.7T 

Z 7 — 1 pi 47T 47T Z 
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where 5, [/, p, and P represent magnetic field, flow speed, density, and pressure, 
the subscript "1" and "2" specify physical quantities in upstream and downstream 
media, and the subscripts "n" and "t" mean the normal components and transverse 
components, respectively. 

The shock solutions of the jump conditions have three possibilities: slow shocks, 
intermediate shocks, and fast shocks, which correspond to three modes of waves in 
MHD. In this thesis we only discuss fast shocks. At a fast-mode shock, the flow is 
decelerated and compressed. The transverse component of magnetic field increases 
as the magnetic field get compressed at the shock. Fast shocks are most frequently 
observed shocks in the heliosphere, including planetary bow shocks, CME-driven 
shocks, most of interplanetary shocks, and the solar wind termination shock. 

Strong shocks in the heliosphere are usually supercritical shocks where the shock 
Mach numbers Ma > M rrtfa >„/ and M,rni-„i < 2.76 dep ends on the shock normal 
angle and the plasma beta, etc ( Stone fc Tsurutanil . ll985l ). In these shocks a fraction 
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Figure 1.5: The structure of quas i-perpendicular supercritical shocks. The picture 
is adapted from (jWu et al.l . Il984j ) with permission by Springer Science + Business 
Media. 



of ions get reflected at the shock front, which provide a n additional diss i patio n 



mechanism as resistivity cannot provide enough dissipation ( iTidman &: Kralll . Il97ll ). 
At small spatial scales, shocks have microscopic structures and the structures are 
often distinct dependent on the shock normal angle Qbu- For quasi-perpendicular 
supercritical shocks, the shock structure is relatively well-defined. Since magnetic 
field vector is mainly perpendicular to the shock normal vector, after ions first 
encounter the shock, their gyro-motions are confined to the vicinity of the shock 
front. The quasi-perpendicular sho cks are featured by the "foot-ramp-overshoot" 
structure, as shown in Figure 11.51 (ILeroy et al.l . |l982j). For quasi-parallel shocks, 
their micro-structures are less clear than that of perpendicular shocks. The micro- 
structure for quasi-parallel supercritical shocks is shown in Figure 11.61 Since the 
background magnetic field is mainly along the shock normal, the reflected ions can 
travel far upstream. It forms ion beams that excite ion-scale low-frequency waves. 
These waves grow in amplitude and is shortened in wavelength as they approach 
the shock. The structures are so called SLAMS (stands for Short Large Amplitude 
Magnetic Structures). 
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Figure 1.6: T he structur e of q uasi-parallel supercritical shocks. The picture is 
adapted from f Treumannl . 120091 ) with permission by Springer Science + Business 
Media. 



1.5 The Motions of Charged Particles 

The heliosphere provides a natural laboratory to study the physics of charged par- 
ticles. In this section we discuss the motions of charged particles influenced by a 
variety of effects. For energetic charged particles in the heliosphere, their motions 
are completely dominated by the Lorentz force. 

When the kinetic energy density of charged particles for the problem of interest 
is much less than the background field, the motion of charged particles has virtually 
no feedback to the background field and can be approximated as test particles. The 
major forces acting on a charged object include electromagnetic force and gravity. 
In spherical coordinates with radial heliocentric distance r, according to Newton's 
second law, we have 

„ dr> 1 GM„m . , 

Ftotai = -± = g(E + -v x B) ^—r + other forces, (1.8) 

at c r z 

where E and B are electric and magnetic field vectors, p, v, q, and m are the mo- 
mentum vector, velocity vector, electric charge, and mass of a particle, respectively, 
c is the speed of light in vacuum, G is the gravitation constant, and M s is the mass 
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of the Sun. 

It is easy to show that for energetic charged particles such as ions and electrons 
in the heliosphere, their motion is dominated by the Lorentz force. For larger parti- 
cles, such as charged dust grains, other effects like radiation pressure and Poynting- 
Robertson drag, etc. may also be important dependent on their charge and mass. 
In this dissertation we focus on energetic charged particles. Their motions in elec- 
tromagnetic field can be described by 

F M = | = gE+^vxB. (1.9) 
at c 

In the simplest case with constant magnetic field B = B z and no electric field, 
the equation has the solution that describes a gyro-motion around a magnetic line 
of force: 

v x — v ± cos(fit + O ) 

v y = — v± sin (fit + (f) ) (1-10) 
v z = v\\ 

x = sin(fit + O ) + x 

v I 

y = — cos(fit + O ) + Vo (l.H) 

Z = V\\t + Zq, 

where v± and vn are velocity components perpendicular and parallel to the magnetic 
field, respectively. The gyrofrequency fi = qB/mc. <po, xo, yo and z are constant. 

When charged particles are moving in electric and magnetic fields that are spa- 
tially and temporally dependent, the motions of energetic charged particles are very 
complicated in general. An important approximation is when the particle moves 
in a slowly varying electric and magnetic field on spatial and temporal scales much 
larger than the scale of gyro-motions. In this case the motion of a charged particle 
can be expressed as a summation of a gyro-motion and a drift motion of the guiding 
center. In a static electric and magnetic field, the guidin g-center motion of charged 



particles can be expressed as (jBoyd fc Sanderson! 120031 ) 
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V 3C 



t; " + ^ B ' (VxB) 



| + Vfi ow x + x VB + x (B ■ V)B, 



(1.12) 

where W± = mv\/2 and W\\ = are the components of particle kinetic energy 

perpendicular and parallel to the magnetic field, and U fi ow ± is the background flow 
speed component that is perpendicular to the magnetic field. The first two terms 
describe the particle motion parallel to the magnetic field including the original 
velocity parallel to the magnetic field and a small modification caused by parallel 
drift. The remaining terms are drift motion transverse to a magnetic field line caused 
by drift in the motional electric field, gradient B drift, and curvature B drift. 

Magnetic fluctuations have profound influences on the behavior of energetic 
charged particles. One important effect is pitch-angle scattering by resonant in- 
teractions between charged particles and magnetic fluctuations. When the resonant 
condition k\\v n — u — = is satisfied, the particle can strongly interact with mag- 
netic fluctuations and the pitch-angles of charged particles may be changed. Here kn 
and u are the wave number parallel to the magnetic field and angular frequency of 
the wave, and v\\ is the parallel velocity of the particle. The resonant wave-particle 
interaction results in the evolution of distribution function that can be described by 
a pitch-angle diffusion. 

1.6 The Transport and Acceleration of Charged Particles 

Particle transport and acceleration are fundamental problems in space physics and 
astrophysics. This section gives an overview of this subject. The large-scale behavior 
of energetic charged particles is usually approximated by a diffusion process, given 
the fact that the scattering time is short compared to the time scale of interest. For 
particles moving in a compression/expansion flow, energetic particles experience an 
increase/decrease in energy. A complete equation that describes the evolution of a 
nearly isotropic dist ribution of en ergetic charged particles is the well-known Parker 



transport equation f lParkerl . Il965l ). The equation describes the large-scale evolution 
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of the distribution function f(xi,p,t) of the energetic particles with momentum p 
dependent on the position Xi and time t including effects of diffusion, convection, 
drift, acceleration and source particles: 



_d_ 

dxi 



dx 



j J 



n di_ v oi.im 

u i o v d,i "Too 

OXj OXi 6 OXi 



df 
dlnp 



Q, 



1.13) 



where is the symmetric part of the diffusion coefficient tensor, Ui is the velocity 
of plasma fluid, and Q is a local source. The drift velocity can be formally derived 
from the drift motion for a single particle in guiding center approximation (Equatio n 
II . 12f) by assuming a quasi-isotropic distribution function ( llsenberg &: Jokipiil . Il979l ). 
It is given by V rf i = (p i cw i /3q i )'V x (B/£> 2 ), where Wi is the velocity of the particles, 
c is the speed of light in vacuum, and qi is the electric charge of the particles. The 
drift can be included in the diffusion tensor as an anti-symmetric part 



k a = p i cw i /?,q i B. 



1.14) 



The motions of energetic charged particles parallel and perpendicular to mag- 
netic field directions are generally quite distinct. The spatial transport coefficient 
along the me a n mag netic field can be related to the pitch- angle diffusion coefficient 
Dfifi (jJokipiil . Il966l ). The diffusion of energetic particles transverse to the mean 
magnetic field is less understood. Recent test-particle simulation has shown that 
the perpendicular diffusion coefficient k i can be a few percent of the parallel diffu- 
sion coefficient Ku ( iGiacalone &: Jokipiil . 119991 1 . The symmetric part of the spatial 
diffusion tensor can be related by the magnetic field vector Bi and diffusion coeffi- 
cient parallel and perpendicular to the magnetic field: 



K±Si 



(k± - K\\)B i B j 



1.15) 



The Parker's transport equation also describes the acceleration (deceleration) of 
charged particles. This equation states that the first-order energy change is due to 



the compr e ssion 



( IKrymsky 



1977 



expansion) of plasma flows. It has been shown by a s e ries p apers 



Axford et al.l . 



1977 



Bell. 



1978 



Blandford fc Ostrikeri . Il978l ) that 
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collisionless shocks are the acceleration sites of charged particles because of the sharp 
compressions. It should also be noted that part icle acceleration may also happen 
in M HD turbulence fe.g.. 



et al 
Wright 



Petrosian fc Liu 



20041). reconnection regions fe.g.. Drake 



2010). and parallel electric fields due to non-ideal MHD effects f Damiano & 



20051 ). These processes will not be discussed in this dissertation. 



1.7 Summary of the Following Chapters 

In this dissertation, we study the acceleration and transport of charged particles, 
focusing on the effect of fluctuating magnetic fields. The importance of turbulent 
magnetic fields has been recognized in many previous works. However, the previous 
studies often treat the effect of magnetic field as a rather simplified process. For 
example, the propagation of energetic particles in space is often assumed to be a 
simple spatial diffusion, and the particle acceleration at collisionless shock has often 
been considered as a process in a 1-D planar shock, etc. While these simplified ap- 
proximations were successful in describing the physical processes, these theories have 
been facing some difficulties in understanding the physical processes and explaining 
observations. 

Recently, there have been a few observations that begin to challenge this picture. 
For e xample, some spacecraft have observed SEP events in great details f Mazur 



et al., l2000l ) . The observations of impulsive SEP events show fine structures in 
intensity-time plots on small temporal sca l es (hours) that are not des cribed by a 



large-scal e spatial diffusion (IMazur et al. 



2000 



Chollet et al. 



2007 



Chollet & 



Giacalone, 120081 ). The Voyager spacecraft crossed the solar wind termination shock 
but did not observe a saturated energy spectrum of ACRs, which is predicted by the 
1-D, time-steady shock theory. As we will show in this dissertation, many features in 
the observations of energetic particles can be explained by considering the turbulent 
nature of magnetic fields. Moreover, we also discuss situations that the transport 
of charged particles cannot be described by the Parker's equation. An example 
is the process related to low-energy particles or other charged particles with high 
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anisotropies. Since the Parker's equation only concerns a quasi-isotropic distribution 
of charged particles, their physical processes can not be appropriately described by 
the equation. Using numerical simulations, we study the transport of energetic 
particles with high anisotropies and the acceleration of low-energy particles. 

In Chapter 2 we study the propagation of charged particles in a turbulent mag- 
netic field, which is similar to the propagation of impulsive SEPs in the inner he- 
liosphere. The trajectories of energetic charged particles in the turbulent magnetic 
field are numerically integrated. The charged particles reached 1 AU are collected 
to mimic spacecraft observations. We show that small-scale variations in the ob- 
served particle intensity (the so-called "dropouts" ) and velocity dispersion observed 
by spacecraft can be well reproduced using this method. Our study also gives a new 
constraint on the error of "onset analysis" , which is a technique commonly used to 
study the propagation of energetic particle and infer the information of the initial 
injection of energetic particles. We also find that the dropouts are rarely produced 
i n the simulations using the so-called "two-component" magnetic turbulence model 
( iMatthaeus et all Il990l ). The result questions the validity of this model in studying 
particle transport. In Chapter 3 we study the acceleration of ions in the existence of 
turbulent magnetic fields. We use 3-D hybrid simulations (kinetic ions and fluid elec- 
tron) to study the acceleration of low-energy particles at parallel shocks. This gives 
new results for the initial acceleration of particles at shocks in fully three-dimensional 
electric and magnetic fields. We also use a stochastic integration method to study 
diffusive shock acceleration in the existence of large-scale magnetic variation. The 
results show that the observations by Voyager spacecraft can be explained by a 2-D 
shock that includes the large-scale magnetic field variation. In Chapter 4 we study 
the electron acceleration at a shock passing into a turbulent magnetic field by using a 
combination of hybrid simulations and test-particle electron simulations. We found 
that the acceleration of electrons is enhanced by including this effect. We discuss 
the application of this process in interplanetary shocks and flare termination shocks. 
We also discuss the implication of this process for SEP events. The correlation of 
electrons and ions in SEP events indicates that perpendicular or quasi-perpendicular 
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shocks play an important role in accelerating charged particles. In Chapter 5 we 
summarize the results and discuss the future work. 
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CHAPTER 2 

The Effect of Turbulent Magnetic Fields on the Propagation of Solar Energetic 

Particles in the Inner Heliosphere 



2.1 Overview of the Transport of Solar Energetic Particles 



One of the most important tasks in the study of solar energetic particles (SEPs) 
is to understand their propagation in the heliospheric magnetic field. The large- 



scale transport of SEPs in the solar wind is usua l 



y stu died by solving the transport 



equation (Equation [LT3]) first derived by iParkerl (119651 ). The spatial diffusion tensor 
(Equation 1 1 . 15|) can be studied 



turbulence (jJokipiil . 



1966 



1971 



jy considering the statistical properties of magneti c 



Giacalone fc Jokipii 



1999 



Matthaeus et al 



2003|) 



The transport equation has been rou tinely us e d to fit the intensity-time profiles of 
SEP events for several decades (e.g., iBurlagal . 119671 ) . For gradual SEP events (see 



Section [L2|) . it has been realized that the profiles of the SEPs cannot be described 
by a simple spatial diffusion process since the energetic particles are continuously 
accelerated at the s hock front, meaning that there must be at least energy changes 
( iKahler et al.l . ll984l ). For impulsive SEP events (see Chapter II .2ft . energetic particles 
are usually released from a confined region in a short duration, which provides a 
relatively simple case to study the transport of energetic particles in space. A long 
standing problem related to the transport of energetic particles is that the mean- 
free paths inferred from SE P events are usually much longe r than those derived 



from the quasi-linear theory (jPalmer 



1982 



Bieber et al. 



19941 ). which assumes that 



the trajectories of charged particles are weakly perturbed by magnetic fluctuations. 
The discrepancy between the observations and the theories is still not well resolved. 
Another problem that requires further investigation is the large-scale transport of 
charged particles normal to the magnetic field. Some analyses give a rather small 
cross-field diffusion so the ratio of the perpendicular diffusion coefficient to the 
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parallel diffusion coefficient is k±/ k\\ ~ 10~ 4 or smaller (IRoelof et all Il983l ). Recent 
numerical simulations and analytical studies give a larger value of k±/k» ~ 0.02 
or larger for energetic part i cles moving in the heliospheric magneti c field at 1 AU 
(IGiacalone fc Jokipiil . 



19991: 



Qin et al 



2002: 



Matthaeus et al. 



2003|). 



Recently, there have been several observations of SEP events that re veal some 
new characteristics of particle transport. For example, Mazur et al.l ( 2000 ) reported 
that the intensity of impulsive SEP events often shows small-scale sharp variations, 
alternatively called the "dropouts" of SEPs (see Figure |2~T|) . These dropouts are 
commonly seen in impulsive SEP events and the typical convected distance between 
the dropouts is about 0.03 AU, similar to the spatial scale of the correlation length 
of the solar wind turbulence. The occurrence of the dropouts does not seem to be 



associated with the rapid magnetic field changes as one can see from Figure 12. 1[ 
meaning that it is more associated with some large-scale transport effects. The 



phenomena indicate that the diffusio n of energetic partic 
magnetic field is very small (see also, 



Chollet fc Giacalone 



es transverse to the local 
and the transport 



2011 



of energetic particles in the solar wind is not a simple diffusion process as described 
by the Parker's transport equation (11.131) . However, Some spacecraft measurements 



indicate that the ratio between perpendi cular to paral 



can reach values of 0.2 or even larger (IDwyer et al 



diffusion coefficien t n i / k \\ 



1997 



Zhang et all 12003 ). 



which is unexpectedly large compared to those obtained from numerical simulations 
(IGiacalone &: Jokipiil . Il999l ). Newly available data shows that impulsive SEP events 
are occasionally seen by all three spacecra ft (STEREO A/B and AC E) with a more 

than 1 00-degree longitudinal separation iwiedenbeck et al. . 2010 ). Giacalone & 

Jokipii ( l2012l )'s numerical simulations suggest that the perpendicular diffusion has to 
be as large as a few percent to explain these multi-spacecraft observations. It should 
be noted that the motions of energetic charged particles transverse to magnetic field 
can be considered to consist of two parts: 1. the actual particle motion across 
the local magnetic field due to drift or scattering; and 2. the particle motions 
along meandering magnetic field lines but normal to the mean magnetic field. The 
observed SEP "dropouts" may be interpreted as that the motions of particles across 
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the local magnetic field is small, but a large part of the perpendicular diffusion can 
be contributed by field-line random walk. These new observations have provided 
an excellent opportunity to examine and constrain the relative contributions from 
these two effects to the large-scale perpendicular transport. 

Using numer ic al sim ulations that contain large-scale turbulent magnetic fields, 
Giacalone et al.l ( 120001 ) have demonstrated that the dropouts can be reproduced 
when energetic particles are released in a small source region near the Sun. The 
idea of the model can be illustrated by Figure 12.21 When the source region of a SEP 
event is small, it just releases energetic particles into a small volume in space so that 
only some magnetic flux tubes are filled by energetic particles. When the field lines 
of force are meandering, the magnetic flux tubes will convect through spacecraft at 1 
AU with a mixture of those filled by energetic particles and those that are not. The 
spacecraft that observe the passage of these flux tubes will see "dropouts" of the SEP 
intensity. It should be noted that this model is consistent with magnetic turbulence 
models that allow a large perpendic ular diffusion with a value of k±/k\\ ~ 0.02 or 
larger due to field line random walk (IGiacalone Sz Jokipiil . Il999l ). The time duration 
between the numerically produced dropouts is several hours, which is similar to that 
observed in the impulsive SEP events. It also naturally reproduces the feature that 
the typical spatial scale for the convected distance betw een the dropouts is the same 
as the correlation scale in the solar wind turbulence flMazur et al.l. 120001). B ased 
on the so-called "two- component" model (see Section l2~2l) . 



Chuychai et al. 



Ruffolo et al. 



((20031) and 



(120071 ) proposed a somewhat different idea. They argued that some 
magnetic field lines in the solar wind can have very restricted random walk. The 
corresponding magnetic flux tubes connecting to the source regions are concentrated 
by energetic particles; For magnetic field lines that are meandering in space, the 
energetic particles in the associated flux tubes will diffuse away. However, this effect 
depends on the "two-component" magnetic field model they use (a composition of 
a two-dimensional fluctuation and a one-dimensional fluctuation) and the motions 
of charged particles during the trapping are not explored by numerical simulations. 
Although previous numerical simulations that contain large-scale field-line random 
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Figure 2.1: An example of imp ulsive SEP events th at show "dropouts" of energetic 
particle intensity (Figure 1 in iMazur et all |2000| ). The figure is reproduced by 
permission of the AAS. 



39 



walk has successfully produced SEP "dropouts" (IGiacalone et al.l . |2000| ). this model 
assumed an ad hoc pitch-angle scattering that is not realistic. Physically, the pitch- 
angle scattering should be caused by s mall-scale magnetic fluctuations, which was 



not present in the 



Giacalone et al. 



(120001 ) model. One main purpose of this study is to 
include the effect of small-scale magnetic turbulence and examine the propagation of 
SEPs in a turbulent magnetic field that has a power spectrum similar to that derived 
from observations. The current numerical simulations directly solve the equations 
of motion for charged particles in turbulent magnetic fields generated by magnetic 
turbulence models. The results (Section [23]) give some new insight for the transport 
of energetic particles in the heliospheric magnetic fields. 

Another important issue on the transport of SEPs in interplanetary space is 
whether we can relate the in-situ spacecraft observation of a particular SEP event at 
1 AU to the initial release of energetic particles in the solar atmosphere (time and/or 
location). Since the energetic particles suffer from spatial diffusion, they gradually 
lose the information about the source regions and injection times after they are 
released. A popular way to get the information about the location of source regions 



and release time is to analy se the onsets of SEP events, name 



particles at a given energy 



2003 



Kahler fc Ragot 



200 



* 



Krucker Sz LinL 



Chollet 



2008 



2000 



Tvlka et al 



Hill et al. 



2009 



y, the earliest arriving 



2003 



Mewaldt et al 



Reames 



2009|). Those 



particles have experienced the least scattering during propagation. One can obtain 
the apparent propagation path length L and the apparent release time t re i ease of 
SEP events by linearly fitting the onsets of the first arriving particles based on the 
formula 



t - t 



release 



L/v, 



(2.i; 



where t is the arriving time for first-arriving particles at a given energy and v 
is the velocity corresponding to the energy. The assumptions implicitly made in 
these studies are that the first- arriving particles are released impulsively and have 
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Figure 2.2: A cartoon that illustrates the effect of meandering magnetic field lines 
on the observations of energetic particles. The SEP dropouts can be seen when 
different field lines pass by the observer. Figure provided courtesy Joe Giacalone, 
University of Arizona. 



41 



experienced no scattering or energy change and that they have travelled exactly 
along the magnetic field lines with pitch-angle cosines \i — 1. A main criticism of 
this method is that the assumption is inconsistent with the fact that the mean- 
free paths of energetic particles in the inner heliosphere are usually less than 1 AU. 
Moreover, the mean- free path is usually energy dependent. Nevertheless, some onset 



analyses do show a good linear relation. In Figure 12.31 we show an example 



or the 



onset analysis to a SEP event, which is adapted from (Figure 8 in iReamesl . 120091 ). 
One can see that the linear relation in the "onset time" versus c/v plot is quite good 
for the energy range they use (> 1 MeV/nucleon). However, the apparent path 
length is larger than typical length of Parker spiral (1.1-1.2 AU). The feature that 
the fitted path length is different from the Parker spiral magnetic field lines has been 



found by a few authors 



Kahler fc Ragot 



2006 



Krucker fc Lin . 



(Met 



2008|). 



2000 



Pei et al. 



ka et al. 



2003 



Mewaldt et al. 



2003 



(120061 ) have demonstrated that the 



effect of large-scale field line meandering can significantly change the arrival times 
for energetic particles, and when some of the field lines are straightened radially, the 
energetic particles can arrive at 1 AU faster than particles travel along the Parker 
spiral. This is the issue that will be further explored in this chapter. Lintunen & 



Vainio (I2004f ) 



Saiz et al 



f l2005h and 



Diaz et al 



(120 111 ) have used more sophisticated 



particle transport models to examine the validity of the velocity dispersion analysis 
and they estimate the errors contained in these analyses could be on the order of 
several minutes or even an hour for typical parameters at 1 AU. However, none of 
these works considers the propagation of energetic particles in a turbulent magnetic 
field that has a power spe ctrum that extends to small resonant scales similar to 



Giacalone fc Jokipiil (119991 ). 



In this study, we use two different types of 3-D magnetic field turbulence models 
often used in studying the transport of energetic particles in space. The generated 
fluctuating magnetic field has a Kolmogorov-like power spectrum with wavelengths 
from just larger than the correlation scale, leading to large-scale field-line random 
walk, down through very small scales that lead to resonant pitch-angle scattering 
of the particles. In Section 12.21 we describe the magnetic turbulence models and 
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Figure 2.3: An example of onset analyses for SEP events. The plot is adapted from 
(Figure 8 in iReamed . I2009I ). The figure is reproduced by permission of the AAS. 
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numerical methods we use to study the propagation of energetic particles. In Section 
12.31 we use the magnetic turbulence models to study the first-arriving particles and 
test the validity of the velocity dispersion analysis. We estimate the errors in this 
technique in a variety of cases with different magnetic fluctuation amplitudes and 
thresholds for the onset. In Section 1274"! we use the magnetic turbulence models to 
study the propagation of SEPs. We show that the "dropouts" of impulsive SEPs 
can be produced using the foot-point random motion model when the source region 
is small. However, for the two-component model, we find that the "dropouts" are 
rarely seen for the parameters we use. The results of this chapter will be summarized 
in Section [2751 



2.2 Numerical Model 

In this study we consider the propagation of energetic particles from a spatially 
compact and instantaneous source in turbulent magnetic fields. We primarily use 
two magnetic field turbulence models that capture the main observations of magnetic 
field fluctuation in the solar wind: the so-called "two- component" model fMat t haeus 



et al.. 119901) and the foot-point random motio n model (e.g 



Jokipii fc Kota 



1989 



Giacalone et al 



Jokipii fc Parker . 



1969 



20061 ). This section gives a mathematical 



description of the turbulent magnetic field models and the numerical method for 
integrating the trajectories of energetic charged particles. 



2.2.1 Turbulent Magnetic Fluctuations 

In a three-dimensional Cartesian geometry (x, y, z), the turbulent magnetic field can 
be expressed as 



B = B + 5B 

= B z + SB x (x,y,z,t)x + SB y (x,y,z,t)y + SB z (x,y,z,t)z. (2.2) 

This expression assumes a globally uniform background magnetic field B and a 
fluctuating magnetic field component <5B. 
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The two- component model is a quasi-static model for the wave-vector spec- 
trum of magnetic fluctu ation based on observations of the solar wind turbulence 
( IMatthaeus et al.l . Il990l ). In this model, the fluctuating magnetic field is expressed 
as the sum of two parts: a slab component 5B S = (B^(z), By(z), 0) and a two- 
dimensional component <5B 2D = (Bl D (x 1 y), By D (x 1 y), 0). The slab component is a 
one-dimensional fluctuating magnetic field with all wave vectors along the direction 
of the uniform magnetic field z, and the two-dimensional component only consists 
of magnetic fluctuation with wave vectors along the transverse direction x and y. 
It has been observed that the magnetic field fluctuation has components with wave 
vectors nearly parallel or perpendicular to the magnetic field, with more wave power 
concentrated in the perpendicular directions (usually about 80% in the solar wind). 
This model captures the anisotropic characteristic of the solar wind turbulence but 
neglects any turbulence component that propagates obliquely to the magnetic field 
B . 

Another often used model for magnetic turbulence is based on the idea that 
magne tic fluctuations can be generated bv foot-point random motions f Jokipii & 



Parker, 



19691 : 



Jokipii fc Kota 



19891 : 



Giacalone et al.l 



20061 1 . One can consider a 



Cartesian geometry with the uniform magnetic field B along the z direction and 
the source surface lying in the x-y plane at z — 0. Since the magnetic field lines are 
frozen in the surface velocity field, the magnetic field fluctuation can be produced 
by foot-point motions in the form of Equation 12.21 We assume that the surface foot- 
point motion is described by Vf p (x,y,t) = V x where \l/ is an arbitrary 
stream function. The fluctuating component of the magnetic field anywhere is given 
by 



SB* = ^ Vfp (x,y,t-z/U). (2.3) 

The magnetic field at z > is assumed to have no dynamical variation but con- 
tinuously dragged outward by a background fluid (the solar wind) with a convective 
speed U. When the magnetic field is evaluated at a certain time, the magnetic field 
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is fully three-dimensional with dependences on x, y, and z. 

In both of these two magnetic fluctuation m odels the magnetic field are variable 



Jones et al. 



i n thre e spatial dimensions. As demonstrated by lJokipii et al.l (119931 ) and 
(119981 ). it is important to consider particle transport in a fully three-dimensional 
magnetic field since particles tie on their original field lines in one- dimensional or 
two-dimensional magnetic field due to the presence of at least one ignorable spatial 
coordinate. The m agnetic fluctuations can be constructed by the random phase 



approximation (e.g., 



Giacalone fc Jokipiil . Il999l ) and assuming a power spectrum of 



magnetic field. This power spectrum can be determined from spacecraft observa- 
tions. The slab component 5B S , two-dimensional component <5B 2£) , and fluctuating 
magnetic field produced by the foot-point random motion 5~B^ P can be expressed as 



5B S = A n [cos a n (cos (p n x + sin (p n y) + i sin a n (— sin (f> n x + cos 4> n y)} 

xexp(ik n z + i/3 n ), (2.4) 



n=l 



<5B 2D = A n i(— sin <p n x + cos 4> n y) 



n=l 



x exp[ifc n (cos (f> n x + sin (j) n y) + i(3 r , 



(2.5) 



Nm / 1 \ 

I j e ik n (cosc/> n x+sm<f> n y)+iu) n (t-z/U)+il3 n ^ 



where (3 n is the phase of each wave mode, A n is its amplitude, u n is its frequency, 
a n is polarization angle, and <p n determines spatial direction of the fc-vector in the 
x-y plane. (3 n , a n , and <p n are random numbers between and 2ir. The frequency 
is taken to be u n = 0.1Uk n . All the forms of fluctuating magnetic field satisfy the 
condition V • 5B = 0. 

The amplitude of magnetic fluctuation at wave number k n is assumed to follow 
a Kolomogorov-like power law: 




Figure 2.4: The turbulent magnetic field lines produced by the foot-point random 
motion model originated from — L c < x < L c and —L c < y < L c at t = 0. See the 
text in Section 12.2.11 for description and parameters. 
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A 2 2 AU n 



II 



(k n L c 



N„ 



n=l - 1 



(2.7) 



where <r 2 = (5B 2 ) jB\ is the total magnetic variance and AV is a normalization fac- 
tor. In one-dimensional, two-dimensional, and three-dimensional omnidirectional 
spectra, AV n = Ak n , 2 r nk n Ak n , and, 4irk 2 Ak n and 7 = 5/3, 8/3, and 11/3, respec- 
tively. 



It has been pointed out by iGiacalone et al.l (120061 ) that these two models are 



closely related and the two-component model can be reproduced using the foot-point 
random motion model by choosing a particular set of fluctuating velocity field. It 
should be noted that both of these two simplified models assume a quasi-static field 
that may not be appropriate for describing magnetic turbulence. Nonlinear struc- 
tures of the magnetic turbulence, such as current sheets that could have important 
effects are not included. Nevertheless, these two models are very useful in studying 
the transport of energetic charged particles in magnetic turbulence and explaining 
the observations of SEP events. We also note that since these models assumes a 
uniform solar wind speed in Cartesian coordinates, they do not include the effects 
of an expanding solar wind in spherical coordinates such as adiabatic cooling and 
adiabatic focusing. 

In the simulations we generally use parameters similar to what is observed in the 
solar wind at 1 AU. The minimum and maximum wavelengths \ m in and \ ma x are 
taken to be 10 _4 -10 -5 AU and 1 AU. The mean magnetic field Bq is typically taken 
to be 5 nT. The convection velocity of the solar wind U is set to be 400 km/s. The 
correlation length is assumed to be L c = 0.01 AU. In figure [2T41 we illustrate 100 
turbulent magnetic field lines origin from a surface region within —L c < x < L c and 
— L c < y < L c at z — at time t — produced by foot-point random motion. It is 
clear that the magnetic field lines are meandering in large scales. The meandering 
field lines originated from the compact region can have large displacements in the x 
and y directions. 
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2.2.2 Test Particle Simulations 

In order to study the propagation of energetic particles in the heliospheric magnetic 
field, we numerically integrate the trajectories of energetic particles in magnetic 
fields generated from the magnetic turbulence models described previously. In each 
time step, the magnetic field is calculated at the position of a charged particle. The 
numerical technique used to integrate the trajectories of energ etic particle s is th e 



so-called Bulirsh-Stoer method, which is described in detail by iPress et al.l ( 119861 ). 
It is highly accurate and conserves energy well. The algorithm uses an adjustable 
time-step method that is based on the evaluation of the local truncation error. The 
time step is increased if the local truncation error is smaller than 10 -6 for several 
consecutive time steps. In the case of no electric field, the energy of a single particle 
in the fluctuating magnetic field is conserved to a high degree with total changes 
smaller than 0.01% during the simulation. 

2.3 A Numerical Study on the Velocity Dispersion of Solar Energetic Particles 

In this section we use both of the magnetic turbulence models described in Sec- 
tion 12.2.11 to study the velocity dispersions of energetic particles in the heliospheric 
magnetic field. The velocity dispersion is due to the fact that faster particles are de- 
tected earlier than slow particles if they are released at the same time and location. 
The parameters are chosen to match the magnetic field observed at 1 AU. Protons 
are released impulsively at z = with random pitch angles \x = v z /v between 
and 1. This injected pitch-angle distribution is different than previous studies f Saiz 



et a!., 120051 ). which assume that the initial pitch angles for all the particles are 
\i — 1. The trajectories of the charged particles are integrated until they reach 
the boundaries at z = 1.2 AU and z = —0.1 AU. The area of the source region is 
taken to be L x x L y = 5L C x 5L C , which is chosen to be larger than the correlation 
length in order to obtain the statistical meaningful results. It is also possible that 
the source regions are small and the transport of SEPs released from those regions 
are only affected by the field lines connecting to the source regions. The effect is 
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unpredictable and requires a demanding computing resource. Here we only discuss 
the situation that the source regions are fairly large. The energies for the released 
protons are 100, 9.4671, 3.3057, 1.6649, and 1 MeV, which correspond to the values 
of 1/v (converted to hour/AU): 0.3, 0.976, 1.65, 2.33, and 3.00, respectively. The 
magnetic variances used in the simulations are varied from a 2 = 0.01 to a 2 = 0.6. 
The mean-fre e paths and parallel diffusion coeff i cients calculated from the quasi- 
linear theory ( Jokipii . 1966 ; Giacalone fc Jokipiil . 1999 ) are listed in Table I2TT1 In 



each case, we numerically simulate the intensity-time profiles for test particles ar- 
rived at 1 AU using the magnetic field generated from four different realizations. 
Each realization is delineated using a new set of random phases, polarizations, and 
propagation angles, etc. The onset times for different thresholds are recorded when 
the values of the intensity reach the thresholds 0.1, 0.01, and 0.001 of the peak 
values, respectively. 

Figure 12.51 illustrates the intensity-time profiles of energetic particles normalized 
using the peak values in the case of the two-component model. In this plot the red 
solid line represents the profile for particles with the energy of 9.4671 MeV, the green 
solid line represents the profile for particles with the energy of 3.3057 MeV, and the 
blue solid line represents the profile for particles with the energy of 1.6649 MeV. 
The thresholds for 0.001, 0.01 and 0.1 of the peak value are labelled using dashed 
lines. Figure 12761 displays a similar plot for the case of the foot-point random motion 
model. The feature of velocity dispersion can be clearly seen from these two plots. 
The propagation of energetic particles along the ^-direction depends on their pitch- 
angles and the scattering they experienced, therefore the energetic particles arrive at 
1 AU at different times. The intensity-time profile usually has a sharp increase when 
the particles start to reach 1 AU. For particles with lower energies, the increases 
are slower compared to the cases for higher energies, presumably because charged 
particles with lower energies have smaller mean-free paths. However, as we will 
show below, for the energy range we simulate (1 - 100 MeV), this usually does not 
introduce a large error in analysing the injection time for energetic particle events 
under the situations that we study. 
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Table 2.1: The par allel mean-free paths and d iffusion coefficients calculated from 
quasi-linear theory (IGiacalone Sz Jokipiil . Il999l ) for particle energies and magnetic 
variances used in this study. 
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Figure 2.5: The intensity-time profiles of energetic particles collected at 1 AU nor- 
malized by their peak values at energies of 9.4671 MeV (the red solid line), 3.3057 
MeV (the green solid line) and 1.6649 MeV (the blue solid line). The magnetic field 
turbulence is generated from the two-component model and the total variance is 



a 



0.3. 
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Figure 2.6: The intensity-time profiles of energetic particles collected at 1 AU nor- 
malized by their peak values at energies of 9.4671 MeV (the red solid line), 3.3057 
MeV (the green solid line) and 1.6649 MeV (the blue solid line). The magnetic field 
turbulence is generated from the foot-point random motion and the total variance 
is a 2 = 0.3. 
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Figure 12.71 displays the intensity-time profiles of energetic particles at 9.4671 
MeV normalized using their peak values in the case of the two-component model. 
In this plot the magnetic variances are a 2 = 0.6 (the blue solid line), a 2 = 0.3 (the 
green solid line), and a 2 = 0.1 (the red solid line). Figure [2781 displays a similar plot 
in the case of the foot-point random motion model. It can be seen that in the case of 
larger magnetic variances, the arriving times for the onsets of energetic particles at 1 
AU are delayed. The delays may be due to two effects: 1. The first-arriving particles 
experience more pitch-angle scattering during the propagation and 2. the lengths 
of turbulent magnetic field lines are larger in the cases of larger magnetic variances. 
Since we study the propagation of particles in a Cartesian geometry, the larger 
magnetic variances only result in longer lengths of magnetic field lines of force. This 
is different than the field lines of force in a spherical geometry, where the Parker's 
spiral field lines can be str aightened and t herefore shortened by the effect of large- 



scale magnetic turbulence (IPei et al. 



20061 ). The decays of the intensities of charged 



particles after the peaks are also different. In the case of larger magnetic variances, 
the decay is slower because of the enhanced scattering during the propagation. 

After obtained the onsets for the intensity-time profiles of energetic particles at 
all selected energies, we linearly fit the onset time t and 1/v based on Equation 12.11 
and get the values of fitted path length L and release time t re i ease . This is similar 
to the method used to analyse the onsets of SEP events. Since we release particles 
at t — 0, the fitted release times actually represent the errors of the onset analyses. 
Two examples of these analyses are given in Figure 12.91 and Figure 12.101 In Figure 
12.91 we present the results of the onset analyses for the two-component model, with 
different magnetic variances varying from a 2 = 0.01 to a 2 = 0.6 and the thresholds 
for intensity onsets are taken to be 0.001 of the peak values. Figure 12.101 shows 
a similar plot for the onset analyses for the cases that the magnetic variance is 
er 2 = 0.6 and the threshold for intensity onset is varied from 0.001 of the peak 
value to 0.1 of the peak value. One can see that the onset analyses usually have 
a good linear relation. For a larger magnetic variance and/or a larger threshold 
value, the onset analyses give relatively large path lengths and more significant 
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Figure 2.7: The intensity-time profiles of energetic particles collected at 1 AU nor- 
malized by their peak values at the energy of 9.4671 MeV. The magnetic field turbu- 
lence is generated from the two-component model and the total variance is a 2 = 0.3 
(the blue solid line), a 2 = 0.1 (the green solid line), and a 2 = 0.03 (the red solid 
line). It is shown that the effect of larger magnetic variance results in delayed onsets 
in the SEP events. 
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Figure 2.8: The intensity-time profiles of energetic particles collected at 1 AU nor- 
malized by their peak values at the energy of 9.4671 MeV. The magnetic field turbu- 
lence is generated from the foot-point random motion model and the total variances 
are cr 2 = 0.3 (the blue solid line), a 2 = 0.1 (the green solid line), and a 2 = 0.03 
(the red solid line). It is shown that the effect of larger magnetic variance results in 
delayed onsets in the SEP events. 
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1.022 


-0.35 


foot-point 


0.03 


0.1 


1.028 


-0.35 


foot-point 


0.1 


0.001 


1.06 


-0.98 


foot-point 


0.1 


0.01 


1.068 


-0.91 


foot-point 


0.1 


0.1 


1.08 


-0.91 


foot-point 


0.3 


0.001 


1.16 


-1.45 


foot-point 


0.3 


0.01 


1.18 


-1.29 


foot-point 


0.3 


0.1 


1.21 


1.29 


foot-point 


0.6 


0.001 


1.26 


-1.29 


foot-point 


0.6 


0.01 


1.296 


-1.48 


foot-point 


0.6 


0.1 


1.34 


-1.09 



Table 2.2: Results of onset time analyses. 
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errors in the release time t re i ease . It has been found that for the case of the foot- 
point random motion model, the estimated errors for onset analyses are smaller than 
those of the two-component model. The reason for this result is probably because the 
parallel diffusion coefficient of particle motion in foot-point random motion model 
is considerably larger than that in the two-component model. This is illustrated in 
Figure 12.171 and will be further discussed in Section 12.41 The results of the onset 
analyses for all cases are listed in Table 12.21 From the table it is shown that the 
errors for the released times t re i ease are usually within several minutes unless the 
variance is large a 2 = 0.6 and threshold = 0.1. This indicates that although the 
pitch-angle scattering could play a role, the onset analyses usually have a small error 
and therefore useful in inferring the release time for energetic particles. However, 
this method is found to have a relative large error in estimating the path length L 
for SEP events. The effect of magnetic turbulence on the apparent path lengths is 
illustrated in Figure [2791 When a larger value of the magnetic variance is chosen, the 
slope of the "l/t> - 1" line is steepened so the apparent path lengths get larger. This 



the typical path l e ngth of the Pa r ker spiral magnetic fie 



Tvlka et al 



20031 : 



Mewaldt et al. 



2003 



Kahler fc Ragot 



1 ( 


Krucker Sz Lin. 


2000 


2006; 


Chollet. 


2008^ It 



should be noted that for a realistic heliospheric magnetic field considering the solar 
rotation, the lengths of the magnetic field lines can occasionally be shorter than 
that of the nominal Parker spiral because some m agnetic fie l d line s are straighten 
radially by the effect of the random magnetic field (IPei et al.l . 120061 1 . It is therefore 



possible that the onset analyses give a path length L smaller than the 



Park er spiral, as thev are seen in some observations ( Hilchenbach et al. 



et al. 



engt 



2003 



is of the 
Chollet 



20071 ). Although we use a different numerical model to test the validity of the 



velocity dispersion analysis , the r esults are qualit atively consistent with previous 



studies ( iLintunen fc Vainio 



2001 



Saiz et al. 



2005|). 
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Figure 2.9: Onset analyses for the two- component model and the threshold is 0.001. 
The different markers represent different magnetic variances a 2 . It is shown that 
the effect of larger magnetic variances change the slope (increases the path length 
L) of the linear fitting in the onset analyses. 
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Onset Analysis (Two component model, o =0.6) 
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Figure 2.10: Onset analyses for the two-component model and the magnetic variance 
is cr 2 = 0.6. The different markers represent different threshold. Their corresponding 
linear fit are labelled using dashed lines. This result shows that effect of the threshold 
on onset analyses. 
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2.4 A Numerical Study of Dropouts in Impulsive SEP Events 

In this section we use turbulent magnetic fields generated from the foot-point ran- 
dom motion model and the two-component mo del to study the S EP "dropouts' 1 
observe d bv spacecraft such as ACE and Wind fjMazur et all 



acalone, 



2008 



2000:Chollet & Gi- 



i). In our test-particle simulations the charged particles are released 
impulsively at z = and their trajectories are numerically integrated until they 
reach the boundaries at z — 1.6 AU and z = —0.1 AU. The spacecraft observations 
at 1 AU are simulated by collecting particles in windows of ci size Lrg x Ly — L c x Lq 
when the particles pass the windows at z = 1 AU. The record for each window is 
plotted as a simulated SEP event observed by spacecraft. The source regions are 
taken to be a circle at the z = plane with a radius of case 1: much smaller than 
the correlation scale (R SO urce = 0.2L C ) and case 2: much larger than the correlation 
scale (R S ource = 5L C ). The energy for the released particles ranges from 20 keV to 
10 MeV. The velocity distribution of released particles is assumed to follow a power 
law / = foV~ A with random pitch angles between and 1. The magnetic variance 
used in the simulations is a 2 = 0.3. In Figure 12.111 we show a simulated SEP event 
using the foot-point random motion model for the case of large source region. The 
upper panel shows the energy-time plot and the lower panel displays the plot of 
inverse velocity 1/v versus the time after the initial release. One can see that in 
this case the simulated SEP event does not show any dropout. In the small source 
region case, the dropout can be frequently seen. An example is given in Figure [2.121 
which illustrates a simulated SEP event energy-time plot (upper panel) and inverse 
velocity 1/v versus time plot (lower panel). It is shown that two SEP dropouts can 
be clearly seen at about t = 13 - 15 hour and t = 17.5 - 20.5 hour, respectively. The 
time intervals of t hese dropouts are ty pically several hours, which is similar to that 



observed in space (IMazur et all 120001 ) 



The velocity dispersion of particles seen by observers can be varied as the par- 
ticles travel along different paths and experience different scattering. This can be 
illustrated by Figure 12.131 which shows an impulsive SEP event plotted as 1/v ver- 
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Time(hour) 



Figure 2.11: An example of SEP event simulated using the foot-point random motion 
model for the case of the large source region. Upper panel: energy-time plot. Lower 
panel: the inverse velocity 1/v versus the time after the release. The simulated 
event does not show SEP "dropouts". 
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Figure 2.12: An example of SEP dropouts simulated using the foot-point random 
motion model for the case of the small source region. Upper panel: energy-time 
plot. Lower panel: the inverse velocity 1/v versus the time after the release. This 
example clearly shows dropouts. 
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sus time. The observation was made by observed by ACE/ULEIS detector in 1999. 
It displays at least two distinct arrival times at 1 AU, which indicates that particles 
follow at least two distinct field-line lengths. In our simulations, we also find that 
in some cases the apparent path lengths can be very different. Two examples are 
presented in Figure 12.141 In these two plots we use blue lines as a reference, which 
represent the particle travel along a field line with length 1.04 AU with pitch angle 
fi = v\\/v = 1. It can be seen from Figure [27141 (upper panel) that the edges of the 
velocity dispersion t = 15 hour and t = 20 hour indicate particles arriving earlier 
than that along the blue line. In Figure 12.141 (lower panel) the earliest arrival time 
for particles at about t = 15 hour and t = 20 are almost the same as indicated by 
the blue line. In addition, the slopes of the edges of the velocity dispersions are 
different. 

We have also attempted to use the two-component model to study the SEP 
dropouts. However, we did not find any clear dropout in our simulations. To further 
resolve this issue, we have prepared two scatter plots that show the positions for 
energetic particles projected on the x - z plane 12 hours after the initial release. 
The results are shown in Figure 12.151 for the foot-point random motion model and 
in Figure 12.161 for the two-component model. It can be clearly seen in Figure 12.151 
that the particles follows the braiding magnetic field lines, and they are therefore 
separated as the field lines doing random walks. However, this feature is not seen 
in Figure 12.161 for the two-component model. A possible reason is that the two- 
component model contain a slab component that can more efficiently scatter the 
energetic particles in pitch-angle. To demonstrate th is, we measure the d i ffusio n 
coefficient by implementing the the technique used by iGiacalone fc Jokipiil ( 119991 ). 
We use the definition of diffusion coefficients = (( 2 )/2t, where ( is the spatial 
displacement in a given time t. We calculate the perpendicular and parallel diffusion 
coefficients for 1-MeV protons in the two turbulence models for the same parameters 
in the simulation. The results are shown in Figure 12.171 For the two-component 
model, we have ku = 1.3 x 10 21 cm 2 /s and k± = 1.1 x 10 19 cm 2 /s. For the foot-point 
random motion model, we have km = 1.5 x 10 22 cm 2 /s and k± = 3.1 x 10 19 cm 2 /s. It 
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Day of Year (1999) 

Figure 2.13: An example of the observed SEP dropouts that show different path 
lengths observed by ACE/ULEIS detector. Figure provided courtesy Joe Mazur, 
Aerospace Corporation. 
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Figure 2.14: Examples of SEP dropouts produced from numerical simulations. The 
result shows a SEP event that observed by two different observers. The two obser- 
vations have different apparent path lengths. 
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is shown that the parallel diffusion coefficient for the two-component model is about 
one order of magnitude smaller than that for the foot-point random motion model, 
meaning that particles experience more scattering in the two-component model. The 
calculation also shows a smaller ratio of k_l/k\\ for the foot-point random motion 
model (acj_//C|| = 0.002) compared with that for the two-component model (k±/k,\\ = 
0.0085). The results question the popularly used "two- component" model. If large- 
scale field line meandering happens to be the explanation for SEP dropouts, pitch- 
angle scattering due to small-scale scattering should be small so energetic particles 
can be mostly confined to their field lines of force. When the pitch-angle scattering is 
large, particles efficiently scatter off their original field lines and observers cannot see 
the intermittent intens ity dropouts. The observat ional evidence of small scattering 
has also been shown by IChollet fc Giacalond (1201 lh . They inferred the intensity-fall- 
off lengths at the edges of the dropouts using ACE/ULEIS data and showed that 
energetic particles rarely scatter off a magnetic field line during the propagation in 
interplanetary turbulence. 



2.5 Summary 

In this chapter, we presented numerical simulations for the propagation of SEPs 
in the inner heliosphere. We numerically integrated the trajectories of energetic 
charged particles in the turbulent magnetic field generated from the commonly used 
magnetic turbulence models. The observations of SEP events were simulated by 
collecting charged particles that reach 1 AU much as a spacecraft detector would. 

Since the initial release of SEPs is highly anisotropic, their propagation cannot 
be well described by the Parker's transport equation. In Section 12.31 we study the 
velocity dispersion of SEPs in the turbulent magnetic field and estimate the error 
involved in the onset analysis commonly used in SEP observations. We find the 
velocity dispersion can be well produced by this model. For a typical turbulence 
variation SB 2 /B^ ~ 0.1 observed at 1 AU and a large source region, we find the 
differences between the apparent release time and the actual release time is less than 
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Figure 2.15: The positions of energetic charged particles projected in x — z plane 
at t — 12 hour. The results are from the numerical simulations using the foot- 
point random motion model. It clearly shows that the particles follows the braiding 
magnetic field lines similar to the Cartoon in Figure 12.21 
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Figure 2.16: The positions of energetic charged particles projected in x — z plane 
at t — 12 hour. The results are from the numerical simulations using the two- 
component model (slab and two-dimensional components). 
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Figure 2.17: Simulated diffusion coefficient of 1 MeV particles parallel (along the 
z-direction) and perpendicular (along the x- and y-direction) to the mean magnetic 
field. Black lines represent results from the two-component model and red dashed 
lines represent results from foot-point random motion model. 
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a few minutes, but the apparent path lengths can be significant different than the 
real path length along the average magnetic field line. For the foot-point random 
motion model, the error for the inferred release time is smaller than that of the 
two-component model. This is due to the parallel diffusion coefficient in the foot- 
point random motion model is considerably larger than that in the two-component 
model. It should be noted that the energetic particles we study on the onset analysis 
is fairly energetic (> 1 MeV) and we assume the energetic particles are released in 
a large source region. In order to increase the statistics, we collect all the particles 
that reach 1 AU rather than collect particles using a window of a certain size. 

We had also reproduced SEP dropouts in the numerical simulations using the 
foot-point random motion model, assuming the SEP source region is smaller than 
the correlation scale. The widths of these dropout are typically several hours, simi- 
lar to the time scales of dropouts observed in space. The velocity dispersion of the 
energetic particles appears to have different path lengths, this indicates that the 
energetic particles travel along different field lines. We have also attempted to use 
the two-component model to numerically simulate the dropouts of energetic parti- 
cles. However, we rarely find the evidence of SEP dropouts in our simulation. This 
is probably because that particle scattering is more efficient in the two-component 
model compare to that in the foot-point random motion model. This result ques- 
tions the popular used "two-component model" in that particles in the turbulence 
model may scatter off field lines too fre quently compared to that constrained by the 
observed dropouts (IMazur et all 120001). This explana tion is also supported by re- 
cent observational analysis ( iChollet fc Giacalonej . 120111 1. which shows that energetic 
particles are mostly confined to their field lines of forces. 

It has been shown that the slab turbulence model gives a me an-free path o f 
energetic particles muc h sm aller than what obs erved in SEP events (jPalmerl . Il982l ). 



Matthaeus et al. 



( 119901 ) and lBieber et al.l ( 119961 ) proposed the two-component model 
yo 2D plus 20% slab) that can give a mean-free path several times larger than 
that given by pure slab model. In this study, we provided the evidence that the 
mean-free path from the two-component model may be need to be even larger to 
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explain the observation of SEP dr opouts 



with the observed mean-free path ( iPalmer 



Mazur et al.1 . |2000| ). which is consistent 



19821 ). 
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CHAPTER 3 

The Effect of Turbulent Magnetic Fields on the Acceleration of Ions at 

Collisionless Shocks 



3.1 Particle Acceleration at Collisionless Shocks: An Overview 



Collisionless shocks in space and other astrophysical environments are considered as 
the main accelerat ors for en e rgetic charged partic l es. D i ffusive shoc k acceleration 



( herein after PSA; iKrvmskv 



Ostriker, 11978 



1977 



Axford et al. 



1977 



Bell . 



1978 



Blandford & 



) is the most accepted theory for the acceleration of charged particles. 



The basic concl usions of PSA can be drawn from the Parker's transport equation 
(Equation II. 13[ iParkerl . Il965l ) by considering a fast shock in a 1-P infinite space and 



steady state system. In the shock frame, we assume that the plasma comes from 
upstream (x < 0) at a flow speed U\, gets compressed and decelerated at the shock 
(x = 0), and flows into downstream (x > 0) at a speed of LAj. For a source function 
Q = $(p ~~ Po)3( x )> which represents the injection of low-energy charged particles 
with momentum po at the shock, the solution of the Parker's equation is 



f(x,p) = (— ) 7 exp( ) x<0,p>p 

PO K xx,l 

(-)~ 7 , x>0 lP >p (3.1) 

Po 

where U\ is the upstream flow speed in the shock frame and k xx .i is the diffusion 
coefficient of the charged particles normal to the shock surface. The solution natu- 
rally predicts a power-law distribution function / oc p' 1 with 7 = 3L r i/(L r i — t/ 2 ). 
For strong shocks with compression ratios between 2.0 and 4.0, the slope index of 
the power law of the distribution function is between 6.0 and 4.0, close to energetic 
particles observed in many different regions of space. Upstream of the shock front, 
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DSA predicts an exponential drop as a fun ction of the distanc e from the shock 



Kennel et al 



1986|). 



similar to some spacecraft observations (e.g., 

Although DSA has been very successful in explaining the acceleration of charged 
particles, this theory does have some difficulties. One of the greatest concerns 
about DSA is how a population of low-energy charged particles gets pre-accelerated 
at collisionless shocks. Since the Parker's equation does not consider low-energy 
particles with high anisotropies, how low-energy particles get accelerated at shocks 
is unclear. This is often referred to as the "injection problem" . There is no current 
consensus on this issue. We will discuss the injection problem in more detail in 
Section 13.21 

Another important problem for DSA is that the observed energetic particles 
associated with shock waves are irregular and variable, and they are sometimes very 
different from what is predicted by the 1-D, steady state solution of DSA (Equation 
13. ip . In contrast, DSA is usually considered to be a simple and robust process. 
It is important to consider how DSA could explain such observations. In recent 
years, it has been realized that the effects of shock geometries, seed particles, and 
spatial and temporal variations, etc. can be important and they are considered to 
be possible solutions for explaining the observed energetic particles. In Section 13.41 
we will discuss these observations and list the possible modifications to the 1-D, 
steady state solution for DSA to interpret the observations. 



There are some other shock-acceleration mechanisms o 



erature. For example, shock drift acc eleration 
and shock surfing acceleration (SSA, 



(SDA, e^ 



Lee et al. 



1996 



ten discussed in the lit 



Armstrong et al. 



Zank et al. 



1985|) 



19961 ) at quasi- 



perpendicular or perpendicular shocks. In SDA, charged particles drift because of 
the gradient of the magnetic field at the shock front. The direction of the drift 
motion is in the same direction as motional electric field vector E = — U x B/c, 
and the particles gain energy during this drift motion. In SSA, it is thought that 
the cross shock potential electric field could reflect ions upstream and the ions gain 
energy by the gyro-motion along the motional electric field. 

SDA and DSA are usually considered to be distinct and their relationship requires 
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some clarification. It has been demonstrated by iJokipiil ( 119821 . 119871 ) that in a 
diffusive process, SDA can be unified into DSA by considering the drift term in 
the Parker's transport equation. In that case, both drift and diffusi on play a role 
and their relative contributions depend on the shock normal angles. IJokipiil (119871 ) 
showed that the acceleration rate is greatly enhanced at perpendicular shock or 
highly oblique shocks since the perpendicular diffusion coefficient is usually much 
smaller than the parallel diffusion coefficient. 

In Section 13.21 we discuss the injection problem of ions for DSA. In Section 13.31 
we present a study on the acceleration of thermal ions at parallel shocks using 3-D 
hybrid simulation (kinetic ions and fluid electron). In Section 13741 we will discuss a 
variety of effects that could modify the simple 1-D diffusive shock acceleration model. 
In Section [3751 we present a study on particle acceleration at shocks containing large- 
scale magnetic variations. 



3.2 The "Injection Problem" (The Acceleration of Low-energy Particles) 



Since the Parker's transport equation (Equation [T7T3], 



Parker 



19651 ) assumes a quasi- 



isotropic distribution of energetic particles, it does not include the pre-acceleration 
process for low-energy particles known as the "injection problem". It is usually 
thought that when the charged particles are energetic enough, they can be efficiently 
scattered by magnetic turbulence close to the shock and get accelerated in DSA. At 
present there is no consensus on this issue. However, one can work out the condition 
for the transport equation to be valid when the anisotropy due to diffusive streaming 
is sm all. The absolute value of the streaming anisotropy is (IGiacalone fc Jokipiil . 
19991 ) 



3\Si\ 3C/i 



vfo 



{K A /Kn) 2 sm 2 9 Bn 



kJk\\ ) 2 sin 2 9 Bn cos 2 9 Bn 



1/2 



.(3.2) 



(ac_l/ac||) sin 2 9 Bn + cos 2 9 Bn 

If we define that the particles can be efficiently accelerated by diffusive shock accel- 
eration when £ < 1, this equation gives an injection velocity 
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v inj > 3?7i < 1 



{k a /k\() 2 sin 2 9 Bn + (1 - K±/K\\) 2 sm 2 9 Bn cos 2 9 m 



[k±/k\\) sin 2 9 Bn + cos 2 9 Bn 



1/2 



(3.3) 



For parallel shocks or quasi-parallel shocks, this indicates that Vi n j > 3U\ and the 
injection is relatively easy. It is usually thought that Alfven waves excited by the 
streaming of high-energy, shock accelerated protons can scatter the pitch angle of 
the particles. When the fluctuations efficiently interact with the particles, these 
particles are trap ped nea r the s hock and gain energy from the plasma compression 



across the shock. lEllisonl (119811 ) first advocated a model for DSA that includes the 
injection process, where the particles are assumed to be originated from the shock- 
heated ions and leak freely from downstream to upstream of the shock. The similar 
mod els have been developed and extended by a number of authors fe.g.. Ellison 



et al, Il990l ; Malkovl . Il998l : Kang et all . |2012[ ). This is usually referred to as the 



However, a number of resear 


chers 


( Quest. 


1988; 


Scholer & Terasawa. 


1990 


; Sc- 


holer, 


1990; 


Kucharek & Scholer. 


1991; 


Giacalone et al. 


1992 


) have found a different 



scenario for the initial energization at parallel shocks based on the results of self- 
consistent hybrid simulations. It is found that the accelerated ions originate from 
the shock layer rather than via leakage from the plasma. The ions can be accel- 
erated from the incident thermal plasma to high energies while they are making 
gyro-motions in the electric and magnetic fields at the shock layer. Although the 
average incident magnetic field is parallel to the shock normal, as the enhanced up- 
stream magnetic fluctuations steepen and convect through the shock layer, the angle 
between the incident magnetic field and the shock-normal right at the shock front 



can be quite large. A particle can gain the first amount of energy by dr i 



being reflected in this "loc ally oblique" shock s t ructu re ( iGiacalone et al. 



ting and 



1992|). It 



has been clearly shown by Kucharek fc Scholer! (119911 ) that most of the accelerated 
part icles are reflected and gain the first amount of energy at the shock layer. Lvu & 



Kan ( 11990l ) have presented the results of hybrid simulations and they claimed that 



7(3 



the leakage protons dominated the accelerated particles. The reason that they ob- 
tained a different result from other researchers may be due to the method they used 
to drive shocks in their simulations. In their simulations, the shocks are initially 
assumed to be a hyperbolic tangent function with a thickness of several ion inertial 
lengths c/upi, where c and u P i are the light speed and proton plasma frequency, 
respectively. The magnetic fluctuations that are important to reflect ions at shock 
front are ignored at the beginning of the simulation. It should be noted that previ- 
ous simulations are restricted to 1-D simulations and occasionally 2-D simulations. 
In those situations the motion s of charged p a rticle s are restricted on their o riginal 



Jones et al. 



( 119981 ). This 



field lines as demonstrated by iJokipii et al.l ( 119931 ) and 
restriction motivates us to study the acceleration of ions at parallel shocks using 
3-D simulations. In Section 13.31 we will present a new study on this problem using 
3-D hybrid simulations. The results show that energetic particles can move across 
field lines but the acceleration mechanism is similar to what is found b y previous 



hybrid si mulations (IQuest 



& Scholer, 



1991 



1988 



Scholer fc Terasawa 



Giacalone et al 



1990 



Scholer 



1990 



Kucharek 



19921 ). Namely, the energetic particles originate 
in the shock layer and are not due to leakage from downstream. 

For perpendicular shocks or quasi-perpendicular shocks, it is thought to be more 
difficult for particles to be injected into P SA. The required pre-ac celeration can be 
achieved by s hock drift acceleration (e.g., lArmstrong et all Il985l ) or shock surfing 



acceleration flLee et al. 



1996 



Zank et al. 



19961 ). Recent progress has been made 



to distinguish the relative importance for these two processes. It has been found 
that in order for shock surfing acceleration to be efficient, the thickness of the shock 



layer has to be very thi n (electron scale) , which is diff erent from w 



rat is observed in 



space (IBale et all 120031 ) and in numerical simulations ( iLeroy et al. 



19821 ). Moreover, 



it has been discussed that the shock thickness has to be fairly large compared to 
electron gyrorad i i to be consistent with the observation of electron heating at shocks 



( iLembege et al. 



20041 ). More detailed st udies using f ull pa rticl e simulations an d 



hybrid simulations have 



jeen presented by I Yang et al.l ( 120091 1 and 



Wu et al. 



(I2009h . 



Giacalone & Jokipii ( 119991 ) has demonstrated that the effect of large-scale mag- 
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netic turbulence can efficiently lower the injection velocity Vi n j at perpendicular 
shocks by increasing the transport of charged particles normal to magnetic field di- 
rection. Recent numerical simulations for the acceleration of charged particles in the 
existence of large-scale magnetic fluctuation s show very efficie nt acceleration, which 
indicates that there is no injection problem ( iGiacalond . 



2005a 



b). It should be noted 



that there is also an injection problem for electrons, which has been considered to be 
more difficult than that of ions. In Chapter H] we will discuss the acceleration of elec- 
trons. The results suggest that the acceleration of electrons prefers perpendicular 
shocks. 

3.3 A 3-D Hybrid Simulation on Particle Acceleration at Parallel Shocks 

In order to examine the acceleration of charged particles at parallel shocks, we per- 
formed 3-D hybrid simulation s to study the initia l acceleration process. As p ointed 



out by a few previous works (IJokipii et al. 



1993 



Giacalone &: Jokipiil . Il994j ) , it is 



important to consider the motions of charged particles in 3-D to avoid the artificial 
restriction on their gyro-motions. 



3.3.1 Numerical Methods 

We perform 3-D hybrid simulations for particle acceleration at parallel shocks. In the 
hybrid simulation (e.g.. lWinske &: Questl . Il988l ). the ions are treated fully kinetically 
and thermal electrons are treated as a massless fluid. It keeps ion-scale kinetic 
physics but neglect electron-scale kinetic physics. Therefore the relevant spatial 
scales are ion gyroradius v t hi/^d and ion inertial length c/u P i, and the relevant time 
scale is the ion gyroperiod fi^ 1 , where v t u is the ion thermal speed, Q ci is the ion 
gyrofrequency, c is the light speed in vacuum, and u P i is the ion plasma frequency. 
This feature is well suited to describe supercritical collisionless shocks, in which the 
dynamics of ions is important. 

In the simulation, the ions are treated as kinetic particles moving in the simula- 
tion domain by solving the equations of motion for each ion j: 
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= qj (E+^^)-r}J (3.4) 

§ - v,, (3.5) 

The first two terms in the right-hand side of Equation 13.41 represent the Lorentz 
force, and the third term describes the effect of resistive coupling between electrons 
and ions where t] is the resistivity and J is the total current. The electric and 
magnetic fields are defined on grid points, and are interpolated to the locations of 
particles. After pushing the particles at each time step, the new particle positions 
Xj and velocities Vj are collected at the grid points to get the quantities such as ion 
density rij, ion flow velocity Vj and ion current Jj = qiUiVi. Since the electrons 
in the simulation are assumed to be a massless neutralizing fluid. We solve the 
momentum equation of electron fluid: 



n e m e -^ = -en e {E + — c ) - V ■ P e + r/J = 0, (3.6) 

where V e is the electron fluid velocity and P^ e is the electron pressure tensor. The 
electron and ion charge density are equal en e = q^i, where n e is the electron number 
density. The pressure tensor P^ in the equation is usually taken to be a scalar 
P^ e = p e T". We assume the electron fluid is adiabatic p e oc n h J 3 . 
The electromagnetic fields are described by Ampere's law, 



V xB = — J=— qMVi-Ve), (3.7) 

c c 



and Faraday's law 



— = -c(V x E). (3.8) 

At each time step, we get V e from Equation 13.71 The electric field is obtained 
from Equation 13.61 and the magnetic field is updated using Equation 13.81 
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We h ave improved the 



models (IGiacalone et al. 



para 



2000 



lelization o 



Giacalone 



;he 1-D, 2-D , and 3-D hybrid simulation 



2004 



2005bl ). The new versions of the 



codes have been implemented and tested on the NASA's Pleiades supercomputer 
using a few thousand CPU cores. In Figure 13.11 we use a 2-D example to illustrate 
the method of parallelization. The 2-D simulation domain is split into several sub- 
domains along the x and y directions. Each sub-domain is placed on a single CPU 
core, CPU 1, CPU 2, ... The sub-domain also has two labels that is determined by 
its location in the global simulation domain. For example, the sub-domain CPU 1 is 
labelled by C x = and C y = 0, and the sub-domain of CPU 6 is labelled by C x = 1 
and C y = 1. The dimension of each sub-domain is d x = L x /n cx and d y = L y /n cy , 
where n cx and Tl C y 3X6 numbers of CPUs along the x and y directions, respectively. 
The information of the kinetic particles in numerical cells Xj and Vj and fields on 
grid points like B, E within each sub-domain is loaded on the corresponding CPU. 
At each time step, the fields at boundaries are transferred to and from its adjacent 
sub-domains (C x + 1, C y ), (C x — 1, C y ), (C x , C y + 1) and (C x , C y — 1). When a 
particle move across a sub-domain, it will be transfer to the corresponding adjacent 
CPU. The number of the CPU is determined by its label C x and C y . 

For this work, we consider a three-dimensional Cartesian grid (x,y,z). All the 
physical vectors have components in three directions and also spatially depend on 
x, y, and z. A shock is produced by using the so-called piston method, in which 
a plasma flow is injected continuously from one end (x = 0, in our case) of the 
simulation box, and reflected elastically at the other end (x = L x ). This boundary 
is also assumed to be a perfectly conducting barrier. The pileup of density and 
magnetic field creates a shock propagating in the — x direction. In the y and z 
direction the boundary conditions are periodic for both particles and electromagnetic 
fields. 

The size of the simulation box for each case is listed in Table 13.11 The Mach 
number of the plasma flow in the simulation frame is Mao — 4.0. The electron and 
ion plasma betas are /3 e = 1.0 and /3j = 0.5, respectively. The grid size is Ax x Ay x 
Az = 0.5c/uj p i x 0.5c/u pi x 0.5c/u pi , where c/u pi is the ion inertial length, and c and 
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Figure 3.1: A 2-D example for the parallization of hybrid simulation. 
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ujpi are the light speed in vacuum and proton plasma frequency in the simulation 
domain at the beginning of the simulation, respectively. The time step is taken to 
be Q C iAt = 0.01, where Q C i is the gyrofrequency. The ratio between light speed 
and upstream Alfven speed is c/v^o = 6000.0, and the anomalous resistivity is r] = 
1 x 10 _6 47ra;~ 1 . The initial spatially uniform thermal ion distribution is represented 
using 25 particles per cell, which is enough to give a reasonable distribution. 

Initially the average magnetic field is assumed to be B = B x, i.e., the averaged 
shock normal angle is 0°. We examine four simulation cases with parameters listed 
in Table 13.11 The spatial sizes in the x direction are taken to be 300c/u P i for all 
the cases. For the first two cases, the system is assumed to have no pre-existing 
magnetic fluctuation but the size of the simulation box in the y and z direction for 
Run 2 (40c/u P i) are considerably larger than that for Run 1 (I0c/u p i). The fluctu- 
ations in the simulation box are self-consistently generated during the simulation. 
In Run 3 and Run 4, the sizes of the simulation boxes are the same as that in Run 
2. Differently, we examine the effect of pre-existing magnetic fluctuations by super- 
posing a random magnetic field on the mean field. The turbulent magnetic field is 
added at the beginning of the simulation and also injected continuously at the x = 
boundary during the simulation. In Run 3, we assume the pre-existing magnetic 
fluctuation is one-dimensional and only depends on x, but the resulting field can 
depend on all three dimensions. The amplitude of the fluctuations at wave number 
k is determined from a Kolmogorov-like power spectrum. The largest and smallest 
wave lengths are taken to be 40c/u p i and 2.5c/ u p i. The total variance is taken to 
be a 2 = 6B 2 /Bq = 0.1. In Run 4 the magnetic power of the injected fluctuation is 
assumed to be isotropically distributed in three wave number vectors k x , k y , and k z . 
The total variance is taken to be a 2 = 0.1, the same as that in Run 3. The difference 
between Run 3 and Run 4 is that in Run 4 the pre-existing magnetic fluctuation is 
allowed to be variable in full three dimensions. As we will see below, this has an 
important effect on the motions of charged particles along the shock front. 
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Run 


L x (c/u pi ) x Ly(c/bj pi ) x L z (c/u pi ) 


o\5B*/Bl) 


Model 


1 


300 x 10 x 10 


0.0 




2 


300 x 40 x 40 


0.0 




3 


300 x 40 x 40 


0.1 


1-D slab 


4 


300 x 40 x 40 


0.1 


3-D isotropic 



Table 3.1: Some parameters for different cases in the 3-D hybrid simulations of par- 
allel shocks: the size of the simulation domain in unit c/u p i, the variance of injected 
magnetic fluctuation a 2 , and the model for the injected magnetic fluctuation. 

3.3.2 Results of Numerical Simulations 

In Figure [3T21 we present the simulation results for Run 1 at VL ci t = 120.0. The figure 
shows the profiles of the magnetic field components in the x direction B x /Bq (a), y 
direction B y /B Q (b), and z direction B z /B (c) along y = 5 c/u C i and z = 5 c/ud, 
where B is the initial upstream magnetic field. It is shown that the magnetic field 
fluctuations around the shock layer are mostly transverse to the initial magnetic 
field in the x direction. The fluctuations are circularly right-handed polarized if 
they propagate along upstream direction, indicating that they are excited by the 
reflected ions that flow upstream. This r esult is consi stent with previous analytical 



theories and numerical simulations (e.g., iQuestl . Il988l ). Close to the shock layer, the 



magnetic fluctuations are steepened and enhanced. The fluctuations get compressed 
and amplified when they are passed by the shock. Different from 1-D simulations, 
the fluctuations in the x direction are also observed. In Figure [3T31 (a) the magnitude 
of magnetic field B/B along y = 5 c/u C i and z = 5 c/u C i is plotted using a black line. 
The magnitude of magnetic field averaged over the y and z directions is overlapped 
on the plot using a blue line. It can be seen that although there is some small 
difference between the profile of the magnetic field and the averaged magnetic field, 
the 1-D cut is very close to the averaged value of the magnetic field, indicating that 
the fields are weakly dependent on the y and z directions. The same features can be 
seen in Figure [3731 (b) and Figure [3731 (c). which show the profiles and averaged values 
of plasma number density h/uq and flow velocity in the x direction V x /Vao, where 
no is the initial upstream plasma number density and Vao is the initial upstream 
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Figure 3.2: 1-D profiles of magnetic field components B x /Bq, B v /Bq, and B z /Bo 
for Run 1 (with a small simulation box in the y and z directions), where Bq is the 
initial upstream magnetic field. 



Alfven speed. 

Figure 13.41 shows the simulation results for Run 2, which has a large simulation 
domain in the y and z directions than that of Run 1. The physical quantities are 
plotted in a way similar to Figure 13.2} except that the profile is selected at y = 35 
c/u C i and z = 35 c/u C i. One can see from Figure [3T41 (a) that the x-component of 
the magnetic field B x /Bq is more fluctuating compared to that for Run 1. This is 
probably due to the fact that a larger number of modes are allowed for the larger 
simulation domain and more obliquely propagating modes are excited. As shown 
in Figure I3.5[ a very different feature for physical quantities in Run 2 compared 
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Figure 3.3: 1-D representation of magnetic-field magnitude, plasma density, and 
velocity for Run 1: (a) the profile of magnitude of magnetic field B / Bq along the 
same line (black line) and the magnetic field magnitude averaged over the y and z 
directions (blue line); (b) the profile of plasma number density n/n along the same 
line (black line) and the number density averaged over the y and z directions (blue 
line), where n is the initial upstream plasma number density; (c) the profile of the 
x-component of the plasma flow speed V x /Vao along the same line (black line) and 
the plasma flow speed in the x direction averaged over the y and z directions (blue 
line), where Vaq is the initial upstream Alfven speed. 
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Figure 3.4: 1-D profiles of magnetic field components B x /B , B y /B , and B z /B 
for Run 2 (with a large simulation box in the y and z directions), where B is the 
initial upstream magnetic field. 
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Figure 3.5: 1-D representation of magnetic-field magnitude, plasma density, and 
velocity for Run 2: (a) the profile of magnitude of magnetic field B/Bq along y = 35 
c/u C i and z = 35 c/w C j (black line) and the magnetic field magnitude averaged over 
the y and z directions (blue line); (b) the profile of plasma number density n/riQ 
along the same line (black line) and the number density averaged over the y and z 
directions (blue line), where n is the initial upstream plasma number density; (c) 
the profile of the x-component of the plasma flow speed V x /Vao along the same line 
(black line) and the plasma flow speed in the x direction averaged over the y and z 
directions (blue line), where Vao is the initial upstream Alfven speed. 
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Figure 3.6: The cross sections of magnetic-field magnitude and density close to the 
shock fronts for Run 1 (x = 130c/u P i) and Run 2 (x = 150c/ 'uj pi ) at Vtdt = 120. 



to that in Run 1 is that the 1-D profiles of physical quantities close to the shock 
front cannot be well represented by the average values over the y and z directions. 
This indicates that the physical quantities in Run 2 is dependent on three spatial 
dimensions, whereas those quantities in Run 1 is only weakly dependent on the y 
and z directions. 

To better illustrate the difference between Run 1 and Run 2, in Figure I37B1 we 
plot the cross sections of magnetic-field magnitude around the shock front for Run 
1 and Run 2. It can be seen that in Run 2 the magnetic- field fluctuation around the 



ss 



shock front has a larger amplitude than that in Run 1. The difference between Run 
1 and Run 2 can also be seen in 2-D contours. In Figure I3~7H we show the color-coded 
contours of the 2-D cuts of (a) the magnitude of the magnetic field B/Bq (the z = 5 
c/u p i plane and the y = 5 c/u P i plane), and (b) the flow velocity in the x-direction 
V x /Vao (the z = 5 c/u P i plane and the y = 5 c/u P i plane) for the case of Run 1. In 
Figure 13.81 we show color-coded contours of 2-D cuts of (a) the density of plasma 
n/riQ (the z = 5 c/u pi plane and the y = 5 c/u pi plane), and (b) the density of 
accelerated particles with energy 3Ei < E < 5Ei (the z = 5 c/u pi plane and the 
y = 5 c/u P i plane) for Run 1, where E\ = \m p Vl is the upstream plasma ram energy 
in the shock frame. The simulation time for these figures is Q C it = 120.0. From these 
figures we can see that the shock front is approximately located at about x = 125 
c/u p i where the physical quantities have a jump. Although the simulation is fully 
3-D, because of the limited size of the simulation box, all the physical quantities 
appear to be close to 1-D and only depend on x. When the size of the simulation 
box is larger in the y and z directions, the restriction on the motion of the particle 
can be removed since the wave modes are allowed to grow in 3-D. This is indicated 
by Figure 13.91 and 13.101 which show the similar contours to Figure 13.71 and 13.81 but 
for Run 2. In this case the lengths of the simulation box in the y direction and the 
z direction are both 40c/cu P j. In Figure [3. Ill the region around the shock front is 
zoomed in to show the magnetic field structure close to the shock. One can see the 
3-D features around the shock front. 

In Figure 13.121 we show a trajectory of a representative particle that is accelerated 
at shock layer for the case of Run 1. The physical quantities for each panel are: (a) 
the energy versus the x position in the simulation frame, (b) the energy versus the 
x position in the shock frame, (c) the energy versus the time in the simulation, (d) 
the time versus the x position in the simulation frame, (e) the time versus the x 
position in the shock frame, and (f) the displacement along the z direction versus 
the displacement along the y direction. One can clearly see that the acceleration 
happens right at the shock front. The acceleration mechanism is due to the reflection 
at the shock layer. As shown in this figure, the particle can ride on the shock front 
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Figure 3.7: 2-D representation of (a) the magnitude of magnetic field B/Bq (the 
z = 5 c/upi plane and the y = 5 c/u P i plane) and (b) the velocity in the x-direction 
V x /Vao (the z = 5 c/u P i plane and the y = 5 c/u P i plane), for the case of Run 1. 
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Figure 3.8: Similar to Figure I3~7T1 2-D representation of (a) the density of the plasma 
flow and (b) the density of the accelerated particles with energy 3E\ < E < 5E\, 
for the case of Run 1. 
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Figure 3.9: 2-D representation of (a) the magnitude of magnetic field B/Bq (the 
z = 20 c/upi plane and the y = 20 c/u P i plane) and (b) the velocity in the x-direction 
Vx/Vao (the z = 20 c/u P i plane and the y = 20 c/w P j plane), for the case of Run 2. 
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Figure 3.10: Similar to Figure 13.91 2-D representation of (a) the density of the 
plasma flow and (b) the density of the accelerated particles with energy 3E\ < E < 
5Ei, for the case of Run 2. 
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and get accelerated for about 15 gyroperiods and the energy gain can be about 
10 times of the plasma ram energy. This process has been discussed by previous 



authors flQuest 



1991 



1988 



Giacalone et al. 



Scholer fc Terasawa. 



1990 



Scholer 



1990 



Kucharek fc Scholer. 



19921 ). Although we use the 3-D hybrid simulation, in this 
case the motion of the particle is still restricted by the limited size of the simulation 
box. This can be seen from Figure 15.121 (f), which shows that the guiding center of 
the particle did not move much during the acceleration process. 

For the case of Run 2, it is expected that the trajectories of the accelerated 
particles could be different since the fluctuations around the shock layer have 3- 
D structures. Figure 13.131 shows an example of the trajectory of the accelerated 
particles for the case of Run 2. We find that the acceleration process is quite similar 
to what is found for the case of Run 1. The main difference is that the particle 
is allowed to cross its original field line in this 3-D electric and magnetic fields. 
This can be seen from Figure [3.131 (f), in which the guiding center of the charged 
particle is drifting in the y-z plane during the acceleration. The difference can be 
seen clearly in Figure 13.141 which shows the three-dimensional trajectories of the 
two representative particles from Run 1 and Run 2, respectively. It can be seen that 
the particle from Run 2 drift along the z direction during the acceleration at the 
shock front, whereas the guiding center of the particle from Run 1 does not move 
much during the acceleration. 

The acceleration of thermal protons at the shock front creates a nonthermal pop- 
ulation of energetic particles. Figure 13.151 shows the density of accelerated particles 
(Upper panel: 3.0E 1 < E < 4.5-E'i; Lower panel: 6.0E 1 < E < lO.OEi) averaged 
over the y and z directions at VL c {t = 120.0, where E\ = l/2rriiUf and U\ is taken 
to be 5.3Vao for both cases. The black solid lines represent the results from Run 1 
and the red dashed lines represent the results from Run 2. In Run 1 more accel- 
erated particles are trapped around the shock front than in Run 2. The energetic 
particles in Run 2 can more efficiently escape upstream and downstream of shock. 
In Figure [3.161 we show downstream energy spectra in 200c/ u P i < x < 300c/u P i at 
VL ci t = 120.0. The black solid line shows the spectra for Run 1 and the red dashed 
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dot line represents the spectra for Run 2. In both of the two fraction of 

thermal protons are accelerated to several times of bulk energy. In the case of Run 2, 
more high-energy protons are concentrated in downstream region and their energy 
range is from several times of bulk energy to about 20 times of bulk energy. 

We also consider the effect of the pre-existing upstream fluctuations. In Run 3 we 
examine the effect of a 1-D pre-existing magnetic fluctuation that only depends on x. 
Figure [3TT7] 13. 18[ and 13. 191 present analyses the same as the previous cases. We find 
that in this case all the physical quantities are close to 1-D. The wave excitation 
processes at the shock front seem to be strongly influenced by the injected 1-D 
fluctuations. Since the injected fluctuation is only dependent on x, it forces the 
motion of the reflected particle to have the one-dimensional dependence, therefore 
the fluctuations excited at the shock front are not strongly dependent on y and z. 
The trajectory analysis (Figure 13.191) of an accelerated particle also confirms this 
idea, which shows that the motion of the particle is close to 1-D and the guiding 
center is approximately localized at the same position in the y — z plane compared 
with that at the beginning of the simulation. 

We further examine the effect of pre-existing magnetic fluctuations in Run 4 
by injecting a 3-D isotropic fluctuation. In Figure 13.201 13.211 and 13.221 we present 
similar analyses to Figure I3.17[ 13.181 and I3.19[ but for the case that the injected 
magnetic fluctuation is isotropic in 3-D. The results show that in this case all the 
physical quantities is dependent on three spatial dimensions. From the analyses of 
the trajectory of a particle, the guiding center drifts in the y-z plane. Note that 
although the guiding center can freely jump to another field line, the acceleration 
process is quite similar to the cases we discussed. 

3.3.3 Summary and Discussion 

The injection problem is a long standing problem for diffusive shock acceleration. 
The charged particles are required to be energetic enough and efficiently interact 
with magnetic fluctuations. In this section we focus on the mechanism for the pre- 
acceleration, in other words, how a population of thermal or low-energy particles get 
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Figure 3.12: The trajectory of a representative particle accelerated at the shock 
layer for Run 1. The physical quantities for each panel are: (a) the energy versus 
the x position in the simulation frame, (b) the energy versus the x position in the 
shock frame, (c) the energy versus the time in the simulation, (d) the time versus 
the x position in the simulation frame, (e) the time versus the x position in the 
shock frame, and (f ) the displacement along the z direction versus the displacement 
along the y direction. 
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Figure 3.13: The trajectory of a representative particle accelerated at the shock 
layer for Run 2. The plotted physical quantities are similar to Figure 13.121 
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Figure 3.14: Three-dimensional plots of representative particle trajectories for Run 
1 and Run 2. In the y and z directions, the displacements y — J v y dt and z = J v z dt 
are used instead of the locations of the charged particles. 
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Figure 3.15: 1-D averages of the density of energetic particles in Run 1 and Run 2 
at n~t = 120.0 





100 




Figure 3.16: Downstream energy spectra for ions in Run 1 (blue) and Run 2 (red) 
at n~t = 120.0 




Figure 3.18: Similar to Figure 13.101 but for the case of Run 3. 
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Figure 3.19: The trajectory of a representative particle accelerated at the shock 
layer for Run 3. The plotted physical quantities are similar to Figure 13.121 
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Figure 3.21: Similar to Figure 13.101 but for the case of Run 4. 



106 




Figure 3.22: The trajectory of a representative particle accelerated at the shock 
layer for Run 4. The plotted physical quantities are similar to Figure 13.121 
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accelerated at the shock front. Previous studies have proposed two different injection 



energization bv downstream heating ( 


Ellison . 


1981 


) an 


i by reflection at the shock 


layer 


Quest. 


1988: 


Scholer & Terasawal. 


199C 


: 


Scholer. 


1990: 


Kucharek & Scholer. 


1991; 


Giacalone et al.. 


1992 


). Although previous self-consistent hybrid simulations 



have found that the initial energization is due to the ion reflection and acceleration at 



the shock layer, the results were ob tained only for 1-D s imulations and occasional 



2-D simulations. As po inted out by 



and 



Jones et al 



Jokipii et al.l (119931 ) . iGiacalone fc Jokipiil ( 119941 ) 



y 



( 119981 ). in a magnetic field that is spatially dependent only on 1-D 
or 2-D, the motions of charged particles are restricted on their original field lines of 
force. In this sense, the previous simulations are not conclusive. 

In this study, we perform 3-D hybrid simulations to investigate the initial accel- 
eration of thermal protons at parallel shocks. In the case of no pre-existing magnetic 
fluctuation, we find that the electric a nd magnetic fields can be generated by the 
plasma instabilities at the shock front (jQuestl . 119881 ). When the size of the simula- 
tion box normal to the initial magnetic field is large enough (40 c/u P i in our case), 
the electric and magnetic fields generated at the shock front show significant 3-D 
variations. By examining the trajectories of the accelerated particles, we find that 
the guiding center of the representative particle can move off its original place when 
the electric and magnetic fields are 3-D, meaning that the particles can move across 
field lines. We also examine if the effect of pre-existing magnetic fluctuations can 
modify this results. We find that for the case of the 1-D pre-existing magnetic fluc- 
tuations, the electric and magnetic fields close to the shock are weakly dependent 
on the y and z directions. Since the injected fluctuation is only dependent on x, it 
restricts the motion of the reflected ions and therefore the fluctuations excited at 
the shock front are not strongly dependent on y and z. When we consider a 3-D 
pre-existing magnetic fluctuation, this restriction is removed. The charged particles 
can move across field lines, as suggested by analysing the their trajectories. In all 
the simulation cases, we find that the charged particles can gain energy at the shock 
layer. The particle can ride on the shock front and gain a large amount of energy 
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(several decades of the plasma ram energy). The results confirm previous hybrid 



simulations that t 



front (I Quest 



1991 



r e initial acceleration of charged particles is right at the shock 



1988 



Scholer fe Terasawa . 



Giacalone et al 



1990 



Scholer 



1990 



Kucharek fc Scholerl . 



19921 ). even if in a 3-D electromagnetic field. 



3.4 Beyond the 1-D Diffusive Shock Acceleration 

Although the 1-D steady state DSA solution gives a very elegant description for the 
acceleration of charged particles at the shock front, some other effects could play a 
role during the acceleration. These effects may help explain the ob s erved variability 
of energetic particles at shocks. For example, Ellison fc Ramaty ( 19851 ) discussed 
the effects such as adiabatic cooling and limited acceleration time, which can cause 
turnovers in the power law energy spectra of shock-accelerated charged particles if 
the diffusion coefficient is proportional to energy. They suggest that the spectra 
have an exponential rollover dJ/dE oc exp(— E/E ). Recent obse rvations indicat e 



that the spectra are more similar to a double power law spectrum ( jMewaldt 



20061 



i.e., af ter the spectrum br eak t he spectra stil l have a power law shape with a steeper 



slope. 



Tylka et al. 



(|2005|) and 



Tylka fc Led (120061 ) argue that the variable spectral 



and compositional characteristics in large SEP events can be produced by consider- 
ing the effects of shock angles and different species of seed particles. Their model 
is based on the argument that the injection speed of particles for diffusive shock 
acceleration is much higher at perpendicular shocks than that at parallel shocks 
(see Equation 13 .3 j) . As a result, parallel shocks can accelerate charged particles 
from solar wind particles, whereas perpendicular shocks only accelerate superther- 
mal particles pre-accelerated in solar flares. This model can qualitatively explain 
the high variable features of spectral characteristics and elemental composition in 
large gradual SEP events. However, recent numerical simulations show that perpen- 
dicular shocks can efficiently accelerate cha rged particles to h igh energy in existence 

2005aJbJ). This questions the 



of strong pre-existing magnetic fluctuation dGiacalone 
validity of the Tylka et al. (2005) model. iLi et al. 



(120091 ) presented a model for 
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particle acceleration at oblique shocks that include both parallel diffusion and per- 
pendicular diffusion. The results can roughly reproduce the observed dependence 
of break energy on the charge-to-mass ratio. It should be noted that most of these 
studies consider particle acceleration at a planar shock. In the real situation the 
particles can sample a range of shock structures with different shock normal angles 
as they move along shock surface. 

Here we mainly discuss a different mechanism that can modify the 1-D solution of 
DSA. DSA is thought to be the mechanism that accelerates anomalous cosmic rays 
(ACRs) in the heliospheric termination shock and also galactic cosmic rays (GCRs) 
with energy up to at least 10 15 eV in supernova blast waves. However, recent in 
situ observations a t the termination shock and in the heliosheath by Voyager 1 
( jStone et al.l . 120051 ) found that the intensity of ACRs was not saturated at the place 
where Voyager 1 crossed the termination shock and kept increasing after entering 
the heliosheath, which strongly indicates that the simple planar shock model is 
inadequate to interpret the acceleration of ACRs. Numerical and analytical studies 
suggest that the possible solution can be made by considering the tempora ry and / or 
spati al variation (Florinski fe Zankl. 120061; McComas Sz Schwadronl. 120061: Jokipii & 



2008; 



Kota fc Jokipiil. 



2008: 



Schwadron et al. 



2008 



). In particular. McComas 



Kota, 

& Schwadron ( 120061 ) discussed the importance of the magnetic geometry of a blunt 
shock on particle acceleration. The idea can be illustrated by Figure 13.231 They 
argued that the missing ACRs at the nose of the heliospheric termination shock is 
due to particle energization occuring primarily back along the flanks of the shock 



where magnetic field 



Kota fc Jokipii 



shock 

resul ts similar to that described bv 



i nes h ave had a longer connection time to the termination 
( 120081 ) present ed a more sophisticated s i mulat ion that gives 



McComas fc Schwadron 



(2006). Schwadron 



et al. ( 120081 ) also developed a 3-D analytic model for particle acceleration in a blunt 



shock, including perpendicular diffusion and drift motion due to large-scale shock 
structure. 

Large-scale magnetic field line meander i ng is ubiquitous in the heliosp h ere and 



other astrophysical environments (1 Jokipii 



1966 



Jokipii fc Parker 



1969 



Parker 
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Figure 3.23: iMcComas fc Schwadronl fl2006l Vs idea on the acceleration of anomalous 
cosmic rays: The Cartoon shows an equatorial cut of the termination shock and 
the heliosphere. The Parker spiral magnetic field get compressed at the termination 
shock. Pickup ions can get accelerated at the termination shock in all directions but 
have different acceleration rate. Because charged particles mainly travel along field 
lines, pickup ions near the nose take a longer time to trav el back to the shock than 
those close to the shock flanks. The figure is adapted from (IMcComas &: Schwadron . 
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Ill 



19791 ). The acceleration of charged-particles in collisionless shocks has been shown 
to be str ongly affected by magnet i c-field turbulence at different scales (IGiacalone . 



2005a 



Giacalone fc Neugebauerl . 120081 ). The large-scale magnetic field variation 



will have important effects on the shock acceleration since the transport of charged 
particles is different in the di rection paralle l and perpendicular to the magnetic 



field, as shown in early work (I Jokipii 



1982 



19871 ). Blunt shocks or shocks with 



fluctuating fronts (ILi fc Zanki . 120061 ) that have a similar geometry, are also relevant 
to this problem. In this study we analyze the effect of the large-scale spatial variation 
of magnetic field on DSA by considering a simple system that captures the basic 
physical ideas. 



3.5 A Study on Particle Acceleration at Shocks Containing Large-scale Magnetic 
Variations 

In this section we present a study on the acceleration of particles at shocks con- 
taining large-scale magnetic variations. The main part of this section has been 
published in the Astrophysical Journal ( iGuo et all l2010l ). Diffusive shock acceler- 
ation at collisionless shocks is thought to be the source of many of the energetic 
particles observed in space. Large-scale spatial variations of the magnetic field has 
been shown to be important in understanding observations. The effects are com- 
plex, so here we consider a simple, illustrative model. Here, we solve numerically 
the Parker transport equation for a shock in the presence of large-scale sinusoidal 
magnetic-field variations. We demonstrate that the familiar planar-shock results can 
be significantly altered as a consequence of large-scale, meandering magnetic lines of 
force. Because the perpendicular diffusion coefficient k± is generally much smaller 
than the parallel diffusion coefficient k», energetic charged particles are trapped 
and preferentially accelerated along the shock front in regions where the connection 
points of magnetic field lines intersecting the shock surface converge, and thus cre- 
ate "hot spots" of accelerated particles. For regions where the connection points 
are separated from each other, the acceleration to high energies will be suppressed. 



112 



Furthermore, the particles diffuse away from the "hot spot" regions and modify the 
spectra of downstream particle distribution. These features are qualitatively similar 
to the recent Voyager's observation in the heliosheath. These results are potentially 
important for particle acceleration at shocks propagating in turbulent magnetized 
plasmas as well as those that contain large-scale nonplanar structures. Examples 
include anomalous cosmic rays accelerated by the solar wind termination shock, 
energetic particles observed in propagating heliospheric shocks, and galactic cosmic 
rays accelerated by supernova blast waves, etc. 



3.5.1 Basic Considerations and Numerical Model 



The diffusive s hock accelerat ion (DSA) can be studied by solving the Parker trans- 
port equation (jParkerl . Il965l ). Here we consider a 2-D system, the Parker transport 
equation (11.131) can be written in Fokker-Planck form as: 
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dp 3 \ dx 

where k xx = (Ax 2 )/2At, k zz = (Az 2 )/2 At, and k X7 = (AxAz )/2At. Following the 
usual approach in stochastic integration (jJokipii &: Levy! . 1 19 771 ). the solution can be 
calculated by successively integrating the trajectories of the pseudo-particles: 



Ax = ri(2K ± At) 1/2 + r 3 (2(«M - K±)At) 
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where r 1; r 2 , and r 3 are different sets of random numbers that satisfy {r^j) = 
and (rf ) = 1. It can be easily demonstrated that the ensemble average of stochastic 
differential equations I3.10[ I3.11[ and 13.121 is the solution of transport equation 13.91 
In order to study diffusive shock acceleration, we approximate the shock layer as a 
sharp variation U x = (U\ + U2)/2 — (U\ — U2)tanh(x/th)/2 with a thickness th much 
smaller compared with the characteristic length of diffusion acceleration k xx i/U\. 
At the same time, we have to make sure that the time step At is small enough to 
resolve the motion in the sho ck laye r. 

Kota & Jokipii (120081 ) and lKotal (120101 ) considered analytically a model in which 
the upstream magnetic field was in a plane (say, x,y), with average direction in 
the y direction. The x-component of the magnetic field was composed of uniform 
sections (straight field lines) alternating in sign, which were periodic in y. They find 
"hot spots" and spectral effects that illustrate the effect of an upstream meandering 
in the magnetic field. 

Here, we consider a 2-D (x, z) system with a planar shock at x = 0, and a 
sinusoid magnetic field B = Bo + sin(fcz)<5B. For most of the parts in this paper 
we discuss the case shown in Figure 13.241 In this figure, the magnetic lines of force 
are illustrated by blue lines. The shock is denoted by the red dashed line. Since, 
in the system of interest here, the magnetic field is small enough that its effects 
are very small (dynamic pressure/magnetic pressure ~ pV^jiB 2 /8tt) ~ 100 in the 
solar wind at 1 AU). The system is periodic in the z direction, with the magnetic 
field convecting from upstream {x < 0) to downstream (x > 0). In the shock 
frame, the particles will be subjected to convection and diffusion due to the flow 
velocities U\ (upstream) and [^(downstream) and diffusion coefficients parallel and 
perpendicular to large scale magnetic field (ftiflin and «2(||,jl))j respectively. The 
gradient and curvature drifts in this case are only in the direction out of the x-z 
plane and thus irrelevant to this study. Because of the steady velocity difference 
between upstream and downstream, charged particles that travel through the shock 
layer will be accelerated. However, since we consider the large-scale magnetic field 
variation, transport of energetic particles in the fluctuating magnetic field become 
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Z 

Figure 3.24: The shock and the magnetic field geometry for an upstream average 
magnetic field perpendicular to the shock normal. The blue lines represent the 
magnetic field lines and the red dashed line indicates the surface of shock wave. 
The flow velocities are U\ (upstream) and U2 (downstream). 



important. The diffusion coefficient in the x-z system can be expressed as: 

Kij = K^bij ^ L . (3.13) 

The normalization units chosen in this study are: the spatial scale X = 10 AU, 
the upstream velocity U\ — 500 km/s, the time scale T = 3 x 10 6 sec, and the dif- 
fusion coefficients are in unit of kq = 7.5 x 10 21 cm 2 /s. The shock compression ratio 
r = U1/U2 is taken to be 4.0. The shock layer is considered to be a sharp variation 
U x = (Ut + U 2 )/2 - (Ui - U 2 )tanh(x/th)/2 with thickness th = 1 x 10" 3 X , which 
is required to be less than k xx i/U\ everywhere in the upstream simulation domain, 
where k xx \ is the upstream diffusion coefficient normal to the shock surface. The 
simulation domain is taken to be [— 2.0X < x < 2.0X , — vrX < z < 7rX ]. The 
parallel diffusion coefficients upstream and downstream are assumed to be the same 
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and taken to be ftp = up = 0.1Xq/T at p — Pq. The ratio between parallel diffu- 
sion coefficient and perpendicular diffusion coefficient is taken to be k±/k\\ = 0.05, 
which is consistent with that determi ned by integrating the traje ctories of test par- 
ticles in magnetic turbulence models fjGiacalone &: Jokipiil . Il999l ). The momentum 
dependence of the diffusion coefficient is taken to be k oc p 4 / 3 , correspondin g to non- 
relativistic particles in a Kolmogorov turbulence spectrum ( iJokipiil . Il97ll ). We use 
the stochastic integration method to obtain the numerical solution of the transport 
equation. In the Equation II. 13[ the source function Q = 5(p — p )5(x) is represented 
by injecting pseudo-particles at the shock x = with initial momentum p = pq. The 
trajectories of pseudo-particles are integrated each time step to obtain the numeri- 
cal solution. The pseudo-particles will be accelerated if they crossing the shock as 
predicted by DSA. The time step is 1 x 10~ 7 T and even reduce at the shock front 
to resolve the variation of U x at the finite shock layer. Particles that move past 
the upstream or downstream boundaries will be removed from the simulation. The 
system is periodic in the z direction, so a pseudo-particle crossing the boundaries in 
the z direction will re-appear at the oppo site boundary and continue to be followed. 



3a) 



2005a) is used in order to im- 



A particle splitting technique similar to (jGiacalond . 
prove the statistics. Although we use very approximate parameters, we note that 
the results are insensitive to the precise numbers. Our results are also qualitatively 
unchanged after allowing the injection rate to vary as a function of shock normal 
angle and different downstream diffusion coefficients. We also note that our model is 
simplified to illustrate the physics of the magnetic-field variation. Other effects such 
as changes in plasma properties (density, temperature, etc.) is small for problems 
of current interest. 



3.5.2 Results of Numerical Simulations 

A shock propagating perpendicular to the average magnetic field 

Consider first the case where the average magnetic field is in the z direction and 
the fluctuating magnetic field is SB = Bq. As shown in Figure 13.241 the magnetic 
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field is convected through the shock front and is compressed in the x direction, 
thus B Z 2 = rB z \. For the sinusoid magnetic field considered in this paper, the 
local angle between upstream magnetic field and shock normal, Obu, will vary along 
the shock surface. As a magnetic field line passes through the shock surface, its 
connection points (the points where the field lines intersect the surface of the shock) 
will be moving apart in the middle of the plane (z = 0) and approaching each 
other on the both sides of the system (|z| = 7rX ). Since, k\\ 3> k±, the particles 
tend to remain on the magnetic field lines. Because the acceleration only occurs 
at the shock front, as the magnetic lines of force convect downstream, the particles 
will be trapped and accelerated at places where the connection points converge 
toward each other, leading to further acceleration. For the regions where the field 
lines separate from each other, the particles are swept away from those regions. 
Figure 13.251 displays the spatial distribution contours of accelerated particles in 
three energy ranges: 3.0 < p/po < 4.0 (top), 8.0 < p/po < 10.0 (middle), and 
15.0 < p/po < 30.0 (bottom). The density is represented by the number of particles 
in simulation and its unit is arbitrary. It can be seen that "hot spots" form in 
the regions that connection points approaching each other at all energy ranges, 
with lobes extend along the magnetic field lines. The density of the accelerated 
particles at the connection-point separating region (in the middle of the plane) is 
clearly much smaller, although there is still a concentration of low-energy accelerated 
particles there since the acceleration of low-energy particles is rapid and efficient at 
perpendicular shocks. At higher energy ranges (middle and bottom), the lack of 
accelerated particles may be interpreted as due to the fact that the acceleration to 
high energies takes time. 

Figure 13.261 illustrates the profiles of the density of accelerated particles for dif- 
ferent energy ranges at z = (top) and z = 7rX (bottom). In each panel, the black 
solid lines show the density of low-energy particles (3.0 < p/po < 4.0), the blue 
dashed lines show the density of intermediate energy particles (8.0 < p/po < 10.0), 
and red dot dashed lines show the density of particles with high energies (15.0 < 
p/po < 30.0). In connection-point separating regions z = (top), it can be seen that 
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Figure 3.25: The representation of density contour of accelerated particles, for en- 
ergy range: 3.0 < p/p < 4.0 (top), 8.0 < p/p < 10.0 (middle), 15.0 < p/p < 30.0 
(bottom). It is shown that the hot spots form on the both side of the system. The 
acceleration at the center of the shock is suppressed. 
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while the downstream distribution of low-energy particles is roughly a constant, the 
density of particles with higher energies increase as a function of distance down- 
stream from the shock. These particles are not accelerated at the shock layer in 
the center of the plane but in the "hot spots". At z — 7rX (bottom), the density 
of particles of all energies decreases as a function of distance, which indicates that 
the accelerated particles diffusive away from the "hot spots" . Since the high energy 
particles have larger diffusion coefficients than the particles with low energy, it is 
easier for them to transport to the middle of the plane. The profile at z = is sim- 



ilar to Voyager's observati on of anoma 



shock and the heliosheath (IStone et al 



ous cosmic rays (ACRs) at the termination 



2005 



Cummings et al.l . 120081 1 that shows the 



intensity of the ACRs is still increasing and the energy spectrum is unfolding over a 
lar ge distance after enter ing the heliosheath. The same physics has been discussed 
by IJokipii fc Kotal (120081 ). where "hot spots" of energetic particles is produced by 
the spatial variation of the injection of the source particles. In our current work, 
the concentration of energetic particles are a consequence of particle accelerated in 
a shock containing large-scale magnetic variation. 

The top panel in Figure 13.271 represents the positions in z direction and the 
times as soon as the particles reached a certain momentum p c = 3po- We show that 
particles are accelerated mainly at the connection-point converging region. There 
are also particles accelerated at middle of the plane because the particles can gain 
energy rapidly at perpendicular or high l y obli que shock due to the smallness of per- 
pendicular diffusion coefficient ( jJokipiil . 119871 ) . However, the further acceleration is 
suppressed by the effect that the charged particles travel away from the connection- 
point separating region (see also the top panel in Figure [3.28H . It is clear that since 
the particles tend to follow the magnetic field lines, when the field line connection 
points separate from each other as field convects through the shock, the particles 
travel mainly along magnetic field and away from the middle of the plane. The char- 
acteristic time for a field line convect from upstream to downstream r c = D/Ui ~To, 
therefore there is no significant acceleration in the middle of the plane after t — r c . 
Some of the particles can get more acceleration traveling from other region to "hot 
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Figure 3.26: The profiles of density of the accelerated particles, for energy ranges: 
3.0 < p/p < 4.0(black lines), 8.0 < p/p < 10.0 (blue dashed lines), 15.0 < p/p < 
30.0 (red dot dashed lines) at different locations z = 0.0 (top) and 7rX (bottom), 
respectively. 
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spots" . Figure 13.271 bottom shows the distance particles traveled in the z direction 
from its original places \z — zq\ versus time when particle get accelerated at a certain 
energy (p c = 3po). It shows that many particles are accelerated close to their original 
position, which is related to the acceleration in the "hot spot" . Nevertheless, there 
are also a number of particles travel from the connection point separating region to 
"hot spot" and get further acceleration, which is represented by the particles that 
travel a large distance in the z direction. 

Figure 13.281 shows the same plot as Figure 13.271 except here the critical momen- 
tum is p c = lO.Opo- It is shown again in Figure 13.281 top that most of particles are 
accelerated to high energy are in the hot spot. However, as opposite to Figure 13.271 
top, there are very few particles accelerated at the center of plane since energetic 
particles more transport away from the middle region and the time available is not 
long enough. For a quick estimate, the acceleration time is approximately, 




(3.14) 



> T c . 



Therefore for most of particles, they do not have sufficient time to be accelerated 
to high energies at the center. A number of particles accelerated at the center will 
travel to the "hot spot" and get more acceleration, as shown in Figure 13.281 bottom. 

Clearly, the acceleration by the shock containing 2-D spatial magnetic field vari- 
ations shows a different picture than the acceleration by a 1-D planar shock, in- 
dicating that the resulting distribution function is spatially dependent. In Figure 
13.291 we show the steady state energy spectra obtained in the regions top: [0.1X < 
x < 0.3X , (tt - 0.2)X < z < 7rX ] (black solid line), and [0.1X < x < 0.3X , 
— O.lXo < z < O.lXo] (green dashed line) and bottom: [O.8X0 < x < l.OXo, 
(tt - 0.2)X < z < ttX ] (black solid line), and [0.8X < x < l.OXo, -0.1X < 
z < 0.1Xq] (green dashed line). It is shown that the spatial difference among distri- 
bution functions at different locations caused by large-scale magnetic field variation 
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Figure 3.27: top: The position in the z direction and time when the particle mo- 
mentum reached p = S.Opo; bottom: The travel distance in the z direction \z — zq\ 
and time when the particle momentum reached p = 3.0po- 
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Figure 3.28: top: The position in the z direction and time when the particle mo- 
mentum reached p = lO.Opo; bottom: The travel distance in the z direction \z — Zq\ 
and time when the particle momentum reached p = lO.Opo- 
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is considerable. The black lines in both top plot and bottom plot, which correspond 
to the "hot spots", show power-law like distributions except at high energies. At 
high energies, the particles will leave the simulation domain before gain enough en- 
ergy that causes the roll over in distribution function, this roll over is mainly caused 
by a finite distance to upstream boundary. For other locations, the 2-D effect we 
discussed will produce the modification in distribution functions. The most pro- 
nounced effect can be found at the nose of the shock (green lines), in the top panel 
the distribution of particles shows a suppression of acceleration at all the energies. 
This insufficient acceleration is most prominent in the range of 6 — 12po- At these 
energies the acceleration time scales are longer than the time for the field line con- 
vection swipe the particles away from the connection-point separating region, as we 
discussed above. The bottom plot shows that deep downstream the spectrum of 
accelerated particles is similar at high energies since the mobility of these particles. 

An oblique shock 

The previous discussion has established the effect of a spatially varying upstream 
magnetic field on the acceleration of fast charged particles at a shock propagating 
normal to the average upstream magnetic field. We next consider the case where 
the shock propagation direction is not normal to the average magnetic field. 

Clearly, if the varying direction of the upstream magnetic field is such that 
at some places the local angle of the magnetic field relative to the average field 
direction exceeds the angle of the average magnetic field to the shock plane, we will 
have situations similar to that discussed in the previous sections. There will be 
places where the connection points of the magnetic field to the shock move further 
apart or closer together. Hence we expect that the same physics can be applicable. 
An example is given in Figure 13.301 In this case the ratio of 5B/Bq is taken to be 
0.5 and the averaged shock normal angle $B n — 70°. It can be seen from this plot 
that the connection points can still move toward each other in some regions. Figure 
13.311 shows the density contours of accelerated particles the same as Figure I3.25[ but 
for the case of the oblique shock. We find that in this case the process we discussed 
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Figure 3.29: The steady state distribution functions at top: [0.1X < x < 0.3X , 
(vr-0.2)X < z < ttX ] (black solid line), [0.1X < x < 0.3X , -0.1X < z < 0.1X ] 
(green dashed line) and bottom: [0.8X < x < 1.0X , (it — 0.2)X < z < rrX ] 
(black solid line), [0.8X < x < 1.0X , -0.1X < z < 0.1X ] (green dashed line), 
respectively. 
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Figure 3.30: The shock and the magnetic field geometry considered for the case of 
an average magnetic field is 70° of the shock normal. The blue lines represent the 
magnetic field lines and red dashed line indicates the position of shock wave. The 
flow velocities are U\ (upstream) and U2 (downstream). 

in the last section is still persistent, even for an oblique shock and relative smaller 
5B/B . The "hot spot" forms correspond to the converging magnetic connection 
points and particle acceleration is suppressed in the region where connection points 
separate from each other. We may conclude that, for a shock that is oblique, if some 
magnetic field lines can intersect the shock multiple times, we have "hot spots" of 
accelerated particles forms where the connection points converging together. 

3.5.3 Summary and Discussion 

The acceleration of charged particles in shock waves is one of the most important 
unsolved problems in space physics and astrophysics. Charged particle transport 
in turbulent magnetic field and acceleration in shock region are two inseparable 
problems. In this section we illustrate the effect of a large-scale sinusoidal magnetic 
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Figure 3.31: The density contour of accelerated particles for the case of oblique 
shock and 5B/B = 0.5, for energy range: 3.0 < p/po < 4.0(top), 8.0 < p/po < 
lO.O(middZe), 15.0 < p/p < 30.0 (bottom). 



127 



field variation. This simple model allows a detailed examination of the physical 
effects. As the magnetic field lines pass through the shock, the connection points 
between field lines on the shock surface will move accordingly. We find that the 
region where connection points approaching each other will trap and preferentially 
accelerate particles to high energies and form "hot sp ots" along the shock su rface, 
somewhat in analogy to the "hot spots" postulated by lJokipii fc Kotal ( 120081 ). The 
shock acceleration will be suppressed at places where the connection points move 
apart each other. Some of the particles injected in those regions will transport to 
the "hot spots" and get further accelerated. The resulting distribution function is 
highly spatially dependent at the energies we studied, which could give a possible 
explanation to the Voyager observation of anomalous cosmic rays. Although we 
have discussed a simplified, illustrative model, the resulting spectra and radial dis- 
tributions show qualitative similarities with the in situ Voyager observations. In 
Figure 13.321 we compare a spatial profile of energetic particle intensity obtained in 
our simulation and the time profile observed by Voyager 1 as it passed the ter- 
mination shock and entered into the heliosheath. One can see that our simulation 
qualitatively agree with the observation. This mechanism gives an interpretation for 
the observation that the ACR intensity did not saturate at the termination shock. 

Thus, the intensities do not in general, peak at the the shock and the energy 
spectra are not power laws. We show that this process is robust even for the case of 
oblique shocks with relatively small magnetic field variations. Large scale magnetic 
field variation, which could be due to magnetic structures like magnetic clouds, or 
the ubiquitous large scale field line random walk, will strongly modify the simple 
planar shock solution. This effect could work in a number of situations for large 
scale shock acceleration including magnetic variations, for example, the solar wind 
termination shock and supernova blast waves. 
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Figure 3.32: A comparison between a spatial profile for the density of energetic 
particles picked from the numerical simulation and observation made by Voyager 1. 
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CHAPTER 4 



The Acceleration of Electrons at Collisionless Shocks 



4.1 Overview 

As we discussed in Section 11.21 and Section 13. 1[ collisionless shocks are efficient ac- 
celerators for a variety of energetic charged particles observed in the heliosphere. 
However, the acceleration of electrons at collisionless shocks is generally considered 
to be more difficult than that of ions. This is primarily due to the fact that the gyro- 
radii of non-relativistic electrons are much smaller compared with that of protons of 



the same energy (by a factor of ym,i/m e ~ 43), therefore low-energy electrons can- 
not resonantly interact with the large-scale magnetic turbulence or ion-scale waves 
close to the shock front. Despite the theoretical difficulties, energetic electrons are 
commonly observed in SEP events. The accelerated electrons are usually observed 
to be associated with quasi-perpendicular sh ocks in interplanetary shocks and plan- 



etary bow shocks (e.g.. lAnderson et al.l . 119791 ). The acceleration of electrons remains 
poorly understood. 

In this chapter we study numerically the acceleration of low-energy electrons 
at shocks that propagate through turbulent magnetized media, using a combina- 
tion of hybrid simulations (kinetic ions, fluid electron) and test-particle electron 
simulations. We find that the acceleration of electrons is greatly enhanced due to 
the effect of large-scale magnetic turbulence, which provides a mechanism for ac- 
celerating electron to high energy and injecting electrons into the diffusive shock 
acceleration (DSA). 

We first review the observations of energetic electrons related to collisionless 
shocks in Section 14.21 Then we give an introduction for the previous theoretical 
works on the acceleration of electrons in Section 14.31 In Section 14.41 we discuss 
the numerical method used in this study. Section 14.51 presents the results on elec- 
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tron acceleration at a shock propagating through a turbulent magnetic field. The 
parameters are similar to observations of electrons associated with interplanetary 
shocks and also the solar wind termination shock. The results can also be used to 
understand the acceleration of electrons at CME-driven shocks. Section H~6l presents 
the numerical results with parameters similar to solar flare region. We discuss and 
conclude the results in Section 14.71 



4.2 Observations of Energetic Electrons Associated with Shock Waves 

Energetic electrons are often observed to be associated with collisionless shocks. At 
the terrestrial bow shock and interplanetary shocks, electrons wi th energy up to 



100 keV are o 



Anderson et al. 



ten observed close to quasi-perpendicular shocks (IFan et al. 



1979 



Gosling et al. 



1989 



Oka et al. 



1964 



20061 ). i.e., the angle between 



the incident m a gnetic field vector and the shock normal #b„ is larger than 45 c 



Anderson et al. 



( 119791 1 showed the ISEE spacecraft measurements of upstream elec- 
trons (> 16 keV) at the Earth's bow shock that originated from a thin region close 
to the p oint of tangency betw een interplanetary magnetic field lines and the shock 



surface. 



Gosling et al.l ( 119891 ) pointed out that the accelerated superthermal elec- 



trons were often seen close to perpendicular shocks or quasi-perpendicular shocks, 
rather than in the vicinity of parallel shocks. The energy spectra of electrons close 
to shocks ap pear to follow a power law with a slope index between —3 to —4. Tsu- 
rutani & Lin (119851 ) reported the observations of energetic electrons associated with 
interplanetary shocks showing "spike-like" flux enhancements for energies > 2 keV. 
The spike events were observed at quasi-perpendicular shocks with Qbu > 70°. Some 
shock crossings had no enhancements of energetic el ectrons that we r e repo rted to be 



Simnett et al 



( 120051 ) have pre- 



associated with low shock speeds and/or small Qsn- 
sented data showing that energetic electrons (~ 50-100 keV) are accelerated close to 
interplanetary shocks. They also showed that some accelerated electrons can escape 
far upstream of a nearly perpendicular interplanetary shock. The clear evidence of 
electron acceleration at interplanetary shocks by DSA is rare, but a recent example 
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discussed by IShimada et al.l (119991 ) shows evidence of the importance of whistler 
waves (a high frequency wave that can resonantly interact with low-energy elec- 
trons) close to a quasi-perpendicular shock. The observation by Voyager 1 at the 
term ination shock showed a spike-like enhancement of energetic electrons (D ecker 



et al., 120051 ). Voyager 2 observed an exponential increase upstream of the termi- 
nation shock and rough ly constant downstr eam in the heliosheath, similar to what 
is predicted from DSA (jpecker et al. . 2008 ). Both of the Voyager spacecraft have 
observed that electrons are accelerated to at least MeV range, indicating that the 
termination shock ca n effi ciently accelerat e elect rons. 

Bale et al. ( Il999l ) and iPulupa fc Bald (120081 ) observed electron foreshocks and 
related Langmuir waves upstream of interplanetary shocks. These events appeared 
to be associated with irregular shock surfaces with spatial scales ~ 2 x 10 5 km. 
They proposed that the complex upstream electron events result from large-scale 
irregularities in the shock surface. The large-scale ripples have been detected using 



2005 



Koval fe Szabo 



multi ple spacecraft observations (INeugebauer fc Giacalone 
20101 ) and the magnitude of radii of curvature is in the range of3xl0 5 to 10 7 km with 
an average of 3. 5 x 10 6 km, sim ilar to the correlation length of the interplanetary 
turbulence (e.g.. IColemanl . Il968l ). 

In solar energetic particles (SEPs) that originate from CME-driven shocks or 
solar flares, electrons are frequently observed to be accelerated to 10 keV - 1 MeV. 



It is important to point out that i n SEP events , the e 



obser ved to have tig 
20091 ) . For example 



it correlations ( 



Posnerl. 



2007 



Oliver 



ectrons and ions are often 



2009 



Haggerty fc Roelol 



Haggerty fc RoeloJ (120091 ) reported that 175-315 keV electrons 



are well associated with 1.8-4.7 MeV protons. ICliverl (120091 ) showed that the elec- 
trons (~ 0.5 MeV) and the protons (> 10 MeV) in large SEP events are strongly 
correlated. The correlation depends weakly on isotope ratio (Fe/O), which indicates 
that the effect of source regions on the electron-proton correlation is not important. 
The tight correlation between electrons and ions has also been found in ground 
base level events (Tylka, private communication), where both electrons and protons 
are accelerated to relativistic energies. These observations indicate that a common 
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mechanism may exist in the process of particle acceleration during SEP events. The 
accelerated electrons are believed to prod uce type II and type H I radio bursts in 
the solar corona and interplanetary space ( Nelson fc Melrose . 19851 ). which provides 
a remote probe for electron acceleration. Herringbone structures featured by fast 
electron speeds are often observed in type II radio bursts, which indi cates that elec 



trons are accelerated to about 20 — 200 k e V at CME-driven shocks ([Roberts . 



Cairns &: Robinson 



1959 



1987 



Cane fc White 



1989|). 



In solar flares, the acceleration of electrons can be observed remotely since the 
accelerated electrons are subjected to several radiation processes. Hard X-ray emis- 
sions produced by superthermal elect r ons a bove magnetic loops have been observed 



(IMasuda et al. 



1994; 



mination shocks (IForbes 



Krucker et al 



1988 



20101) and may be associated with flare ter- 



Shibata et al 



19951 ). U sing the em i ssions from 

(boOJ) have 



Shih et al. 



neutron-capture 7-ray line and electron bremsstrahlung, 
shown that the emissions of > 30 MeV protons and that of > 0.3 MeV electrons 
are proportional to each other. Again, this indicates that the electrons and ions are 
accelerated by closely related mechanisms. 

In astrophysical shocks, electrons are observed to be accelerated to very high 
energies (may reach 10 15 eV) and produce stro ng synchrotron emissions in the exis- 
tence of magnetic fields ([Reynolds et al.l . l201ll ). At those energies, the electrons are 
ultra-relativistic and their gyroradii are virtually the same as protons of the same 
energy. In that situation, both the acceleration of electrons and that of protons may 
be considered to be DSA. 

As a summary, there is plenty of evidence of energetic electrons associated with 
fast shock waves. Observations suggest that accelerated electrons are more asso- 
ciated with oblique shocks and rarely seen at parallel shocks. Compared to ions, 
energetic electrons are more-easily to excite electromagnetic radiation (radio bursts, 
bremsstrahlung emission, and synchrotron radiations, etc.), which can be used as 
remote tracers for energetic electrons. In many situations, electron acceleration and 
ion acceleration are observed to be correlated, which may give a constraint to the 
acceleration mechanism. 
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4.3 Review of Previous Theoretical Works 



In or der to explain the energization of electrons within the shock layer. Holman & 



Pesses (119831); 



Wul (119841 ) and iLeroy &: Mangeneyl (119841 ) developed analytic models 



for electron acceleration from thermal energies by adiabatic reflection off a quasi- 
perpendicular shock. This is usually referred to as shock drift acceleration (SDA) 
or fast-Fermi acceleration. The theory describes a scatter-free electron acceleration 
process in a planar, time-steady shock. It obtains a qualitative agreement with 



observations at Earth's bow shock in terms of the loss-cone pitch-an g! 



e distribu- 



tion and energy range of accelerated electrons. iKrauss-Varban et all ( 119891 ) used 
the combination of electron test-particle simulation and 1-D hybrid simulation and 
verified Wu's basic conclusions. The main energy source of fast Fermi acceleration 
comes from the —V x B /c electric field that is the same as shock drift acceleration 
( [Armstrong et all 119851 ). It can be demonstrated that fast-Fermi acceleration and 



SD A are the same process in two different frames of reference f Krauss-Varban & 



Wu, 



19891 ). Thus one would expect electrons to drift in the direction perpendicular 



to the flow and magnetic field during the acceleration at the shock front. 

However it is commonly known that in this process bo th the fraction and at 



tamable energy of the accelerated particles are limited (e.g.. lBall fc Melrosd . l200ll ). 
This cannot explain t he observed high-ene rgy electron. For example, observations 
at Earth's bow shocks (IGosling et all Il989l ) suggested that the accelerated electrons 



have a power law distribution with an expo netial drop, which ca nnot be produced 



by this simple mechanism. In a recent work, iPulupa et al 



(120121 ) presented a com- 



parison between the STEREO measurements of electron acceleration at Earth's bow 
shock and the theory of adiabatic fast-Fermi acceleration (jWul . Il984l ). It is clear that 
the fast-Fermi acceleration fails to predict the observed energetic electrons in both 
the accelerated fraction and spectrum sh ape. Herringbone structures observed in 
Type II burst (jCairns Sz Robinsonl . 119871 ) require a high electron energy, and the 
energetic electrons cannot be produced by fast Fermi acceleration unless electrons 
are extremely hot. 
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Some recent progress in shock acceleration is the consideration of n onplanar 



effects suc h as shock ripples a nd magnetic turbulence. The simulations by [Burgess 



fl2006h and 



Umeda et al. 



( 120091 ) show that small-scale shock rip ples can be important 



Savoini et al 



(2010) 



in scattering the electrons and facilitating the acceleration 
discus sed the non-adiabatic motion of electrons. Recently, iJokipii &: Giacalone 
(120071 ) proposed a novel mechanism to solve the injection problem that does not 
require strong pitch-angle scattering from small-scale fluctuations. The low-energy 
particles can move along the meandering field lines of force, travel back and forth 
across a shock front, therefore gain energy from the difference between upstream and 
downstream flow velocities. In this study (Section 14.51 and Section I4T61) . using self- 
consistent hybrid simulations combined with test-particle simulation for electrons, 
we have found that efficient electron acceleration can happen after considering large- 
scale pre-existing upstream magnetic turbulence. The turbulent magnetic field leads 
to field-line meandering that allows the electrons to get accelerated at a shock front 
multiple times. The shock surface becomes irregular on a yariety of spatia l scale s 
from small-scale ripples due to ion-scale plasma instabilities ( Lowe fc BurgesJ . 2003 ). 



to large-sca le structures caused by th e interaction between the shock and upstream 



turbulence (IGiacalone fc Neugebauer 



2008 



Lu et al 



20091 ). The rippled surface of 



the shock front also contributes to the acceleration by mirroring electrons between 



the ripples. The observa tional evidence of these large-sca 



& Bale, 



2008 



by a n umber of authors (INeugebauer fc Giacalond. 



Koval &; Szabo . 



on acceleration of ions ( 



2005 



e ripples has bee n shown 



Bale et al. 



1999 



.Pulupa 



Giacalone 



201o[). Th ese results, along with the previous work 
2005al ibl). suggest that large-scale turbulence has 



an important effect on the acceleration of both electron and ions at shocks, which is 
consistent w ith the correlation betwee n ions and ele ctrons in solar energetic particle 



events (e.g., 



Haggerty fc Roeloi 



2009 



Oiyerj, I2OO9J) 



Another scenario that solves t 
excited close to the shock front. 



lis inject i on problem relies on small-scale waves 



Levinsonl (11992 



1994J ) has proposed an analyti- 



cal theory for thermal electron injection due to whistler waves excited by electrons 
streaming upstream. In this theory efficient electron acceleration requires large 
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Mach numbers Ma ~ 43/ where f3 e is the ratio of thermal electron pressure to 
magnetic pressu re inside the shock. This condition is rarely satisfied for shocks in 
the heliosphere. lAmano fc Hoshinol (120071 ) proposed "electron surfing acceleration" 
where electron can reach efficient acceleration by being trapped in electrostatic soli- 
tary waves. This mechanism also requires a relat ive high mach number shock (A lfven 
mach number Ma is greater than 14). Recently, iRiquelme fc Spitkovskyl ( 1201 ll ) used 
full particle simulations to study electron acceleration at oblique shocks. For the 
parameters they use, they find that efficient nonthermal electron acceleration when 
electrons are scattered by oblique whistler waves at a Ma ~ 7 quasi-perpendicular 
shock. However, these full particle simulations generally use unrealistic mass ratio 
mi/m e and the ratio between the speed of light and thermal speed c/V t h e , which may 
produce some artificial effects. In fact in the simulations they find the electrons are 
heated at t he shock laver to a temperature the same as that of protons ( Riauelme & 
Spitkovsky, 120111 ) . This is not consistent wi th in situ observ a tions at interplanetary 



shocks and the relevant theoretical works (IThomsen et al. 



1988 



1987 



Schwartz et al. 



Scudderl . Il995l ) , in which the the 



less than that of ions (IThomsen et al 



leafing of electrons at shocks is significant 



198' 



Schwartz et al 



198 



Scudder 



19951 ). 



y 



The effect of whistler waves for electron acceleration at shocks still remains to be 
clarified. 

4.4 Numerical Method 



Investigating the motion of charged particles in the vicinity of a collisionless shock 
requires a spatial scale large enough for particles to travel back and forth across 
the shock, and a spatial resolution small enough to include the detailed physics for 
particle scattering and shock microstructure. We implement a combination of a 2-D 
hybrid simulation to model the fields and plasma flow and a test particle simulation 
to follow the orbits of a large number of ene rgetic elect r ons. I n the first step, we 
employ a two-dimensional hybrid simulation ( Giacalonel . 2005bl ) that includes pre 



existing large-scale turbulence. In the hybrid simulation (e.g., IWinske fc Quest 
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19881 ). the ions are treated fully kinetically and the thermal (i.e., non-energetic) 



electrons are treated as a massless fluid. This approach is well suited to resolve 
ion-scale plasma physics that is critical to describe supercritical collisionless shocks. 
We consider a two-dimensional Cartesian grid in the x — z plane. All the physical 
vector quantities have components in three directions, but depend spatially only on 
these two vari ables. A shock is produ ced by using the so-called piston method (for a 



discussion, see IJones fc Ellison 



199ll ). in which the plasma is injected continuously 



from one end [x = 0, in our case) of the simulation box, and reflected elastically at 
the other end (x = L x ). This boundary is also assumed to be a perfectly conducting 
barrier. The pileup of density and magnetic field creates a shock propagating in the 
—x direction. To include the effect of large-scale magnetic fluctuations, a random 
magnetic field is superposed on a mean field at the beginning of the simulation 
and is also injected continuously at the x = boundary during the simulation. The 
simplified one- dimensional fluctuations have the form B(z, t) = 5~B(z, t) +Bi, where 
B x is the averaged upstream magnetic field. The fluctuating component contains 
an equal mixture of right- and left-hand circularly polarized, forward and backward 
parallel-propagating plane Alfven waves. The amplitude of the fluctuations at wave 
number k is determined from a Kolmogorov-like power spectrum: 



P(k) oc 



1 



(4.1) 



l + (LL c ) 5 /3' 

in which L c is the coherence scale of the fluctuations. For the simulations presented 
in this study, we take L c = L z , which is the size of simulation box in the ^-direction. 
Note that in addition to magnetic fluctuations, there are also velocity perturba- 
tions with 5v = vaiSB/Bi (Alfven waves). Different from previous studies, the 
consideration of large-scale magnetic fluctuations enables us to consider the effect 
of pre-existing magnetic turbulence on electron a cceleration , which has been shown 



to be important for low-energy ion acceleration ( iGiacalond . l2005al Jbl) since particle 
transport normal to the mean field is enhanced. However, particle transport in full 
3-D turbulence cannot be properly treated in a self-consistent way using available 
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compu tation 



Jokipii, 



. As demonstrated by previous work (IJokipii et al 



1994 



Jones et al 



1998 



1993 



Giacalone & 



), in the model with at least one ignorable coordi- 
nate, the center of gyration of particles is confined to within one gyroradius of the 
original magnetic field line. The test-electrons can still move normal to the mean 
field in our model because of the field-line random walk. 

In the second part of our calculation we integrate the full motion equation of an 
ensemble of test-particle electrons in the electric and magnetic fields obtained in the 
hybrid simulations (see Figure |4~T|) . This part of the calculation is done separately 
from the ma in hybrid simulation as a post processing phase. As noted bv Krauss- 
Varban et al. ( 119891 ). high-order interpolation of fields is required to ensure numerical 
accuracy and avoid artificial scattering in calculating electron trajectories. In this 
work we use a second-order spatial interpolation and a linear temporal interpolation, 
which ensure the smooth variations of the electromagnetic fields. The test-particle 
electrons are released uniformly upstream when the shock has fully formed and 
is far from the boundaries. The numerical technique used to integrate electron 
trajectories is the so-called Bulirsh-Stoer method, which is described in detail by 



Press et al. 



( 119861 ). It is highly accurate and conserves energy well. It is fast when 
fields are smooth compared with the electron gyroradius. The algorithm uses an 
adjustable time-step method based on the evaluation of the local truncation error. 
The time step is allowed to vary between 5 x 10~ 4 and O.lfi" 1 , where Q ce is the 
electron gyrofrequency. The ratio Q ce /Q C i is taken to be the realistic value 1836. 
The total number of electrons in the simulation is 1.6 x 10 6 . The electrons that 
reach the left or right boundary are assumed to escape from the shock region and 
are removed from the simulation. The bound ary cond i tion i n the z direction is taken 
to be periodic. The readers are referred to f BurgesJ . 2006 ) for more details on the 
numerical methods. 

Magnetic field turbulence has already proved to have key effects on the particle 
acceleration in collisionless shocks. Unfortunately, solving the whole problem in 
three-dimensional space and resolving magnetic turbulence from coherence scale to 
electron scale are still limited by available computation in the near future. This 
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limitation motivates us to solve these problems approximately. We also note that in 
our model the electron test-particle simulation is not self-consistent since the hybrid 
simulation does not include the electron scale physics. The electron scale shock 
structure may be important but is neglected here. 

4.5 The Effect of Large-Scale Magnetic Turbulence on the Acceleration of Elec- 
trons: Interplanetary Shocks 

Much of this section has been published in the Astrophvsical Journal f Guo & Gi- 



acalone, 



201(1 ) . 



4.5.1 Initial Condition and Parameters 

In this section we use the numerical method described in Section 14.41 to study the 
acceleration of electrons at a shock that propagates through a turbulent magnetic 
field. The parameters are similar to interplanetary shocks and the solar wind ter- 
mination shock. For most part of this study, we consider a turbulence variance 
cr 2 = 5B 2 / B\ = 5v 2 /v A1 = 0.3, where 5v and v^i are the magnitude of velocity 
perturbation and upstream Alfven speed, respectively. We also discuss the effect of 
different values of turbulence variances. The size of the simulation box for most of 
situations is L x x L z = 400c/o;pj x 1024c/u P i, where c/u P i is the ion inertial length. 
We also examine the effect of different size of simulation box. The Mach number of 
the flow in the simulation frame is Mao = 4.0 and the averaged Mach number in the 
shock frame is about 5.6. Most of the results presented here are for the averaged 
shock normal angle (#B n ) = 90°, but we also simulate the cases for (#B n ) = 60° 
and 75° to examine the dependence of the acceleration efficiency on shock nor- 
mal angle. The other important simulation parameters include electron and ion 
plasma beta f3 e = 0.5 and = 0.5, respectively, grid sizes Ax = Az = 0.5c/ 'ui pi , 
time step At = O.Olfi" 1 , the ratio between light speed and upstream Alfven speed 
c/vai = 8696.0, and the anomalous resistivity is taken to be r] = 1 x 10~ 5 47roj~ 1 . The 
initial spatially uniform thermal ion distribution was generated using 40 particles 
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per cell. Under these parameters, the upstream Alfven speed is about 34.5 km/s 
and the shock speed is about 193 km/s in upstream frame. 

We assume non-relativistic motion that is reasonable because the highest energy 
electrons obtained in our study are still non-relativistic. We release a shell distribu- 
tion of electrons with energy of 100 eV, which corresponds to an electron velocity 
V e = 30.7f/i = 5.7v the in the upstream frame, where U\ is upstream bulk velocity 
in the shock frame and v t he is the thermal velocity of fluid electrons considered in 
the hybrid simulations. This energy is typical for the halo component of electron 
velocity distributions observed in the solar wind. 



4.5.2 Simulation Results for the Acceleration of Electrons 



Figure 14.11 shows a snapshot of the z component of the magnetic field, B z /Bi, at 
time 110O" 1 in a color-coded contour. At this time, the shock is fully developed. In 
this case, the angle between the average magnetic field direction and shock normal, 
(9btl) is 90°. The position of the shock front is clearly seen from the boundary 
of the magnetic field jump. The shock is moving in the — x direction at a speed 
dependent on z, which is about 1.6 v&\ on average. On the right bottom a small 
region of the simulation domain is zoomed in, which shows small-scale irregularities 
at shock front. Because of the effect of large-scale turbulence with the shock, the 
shock surface becomes irregular on a variety of spati al scales from small-sca le ripples, 
which could be due to ion-scale plasma instabilities ( iLowe fc Burgesd . 120031 ) . to large- 
scale structure caused by the inte r action between the shock and up s tream turbulenc e 



(INeugebauer fc Giacalone 



2005 



Giacalone fc Neugebauer 



2008 



Lu et al. 



2009|). 



Locally, the structure of the shock, shown in top right panel, is still clearly a quasi- 
perpendicular shock. The upstream magnetic field is compressed and distorted as 
it passes through the shock into the downstream region. We note that the rippling 
of the shock and varying upstream magnetic field leads to a varying local shock 
normal angle along the shock front. As we will discuss later, the irregular shock 
surface and magnetic field geometry will efficiently accelerate electrons and produce 
a number of features similar to observations, such as the electron foreshock and 
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spike-like intensity increases at the shock front. The meandering of field lines close 
to the shock surface helps to trap the electrons at the shock, leading to efficient 
acceleration. The shock ripples also contribute to the acceleration by mirroring 
electrons between them. 

Figure H21 shows a color-coded representation of the number of energetic electrons 
with energies higher than 10 times (i.e., 1 keV) the initial (at release) energy at 
three different times (a) 76fi~ , (b) Slfi" 1 , and (c) 90Q~^. It is found that after the 
initial release, a fraction of the electrons are reflected and accelerated at the shock 
front, and then travel upstream along the turbulent magnetic field lines. These 
accelerated electrons are then taken back to the shock by the field line meandering, 
which provides even further acceleration. The number of energetic electrons close to 
the shock surface is highly irregular because the acceleratio n efficienc y varies along 
the shock front depending on the local shock normal ang le Jwilligsi l. Most of the 
electrons concentrate near the shock front since the global magnetic field is mostly 
perpendicular to the shock normal. As the field lines convect through the shock, the 
electrons eventually are taken downstream. Since the electrons are tied to individual 
field lines in 2-D magnetic field, once the electrons are no longer capable of crossing 
the shock, there will be no additional significant acceleration. At this point, once 
all electrons are downstream, the energy spectrum no longer changes with time. 

Examination of the trajectories of some electrons shows that the rippling of 
the shock front also contributes to the acceleration by mirroring electrons between 
the ripples, as illustrated in Figure 14.31 In this figure, the top left plot displays 
the trajectory of a representative electron in the x-z plane, overlapped with the 2-D 
gray-scale representation of B z at Vt c {t = 89.0, the gray scale is the same as in Figure 
14.11 The upper right plot shows the position of this electron (in x) as a function of 
time. The electron bounces back and forth between the ripples for several times. For 
example, the reflections are labeled a-b, c, d, f-g, h, and j. The energy change as 
a function of position, x, corresponding to these reflections is shown in the bottom 
left panel. We find that there are jumps in energy at each of the reflections. The 
panel on the bottom right shows the electron energy as a function of time that also 
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Figure 4.1: A snapshot of the magnetic field in the ^-direction B z jB\ represented 
in color-coded scale at t — HOf^ 1 , where B\ is the averaged upstream magnetic 
field strength. A region at the shock front is zoomed in on the right bottom and 
a profile is illustrated on the right upper panel. The shock surface is shown to be 
rippled and irregular in different scales. 
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Figure 4.2: The number of energetic electrons with energies > IOEq, the initial 
release energy E = 100 eV, at (a) fl ci t = 76, (b) fldt = 81, and (c) fl C it = 90, 
respectively. Initial electrons are released uniformly upstream at Q ci = 70. 



illustrates the features of multiple accelerations related to multiple reflections. The 
trajectory analysis shows that the electron will be mirrored between the ripples and 
get accelerated multiple times. Note that the shock does not move much during the 
time scale of this trajectory. 

We now consider the effect of varying the angle between the mean magnetic field 
and shock-normal. Shown in Figure 14.41 are the resulting energy spectra for three 
different mean shock-normal angles ({0Bn) = 60°, 75°, and 90°, respectively) at the 
end of simulations (Q C it = 120.0). It is found that for (0Bn) = 90°, the electrons 
can be readily accelerated to up to 200 — 300 times the initial energy within SOfi" 1 . 
The spectrum is flat between abo ut 0.1 keV to .7 keV. This shape is similar to the 
"plateau" structure discussed by iBurgessl (120061 ). Above 1 keV, the spectrum falls 
off with energy with a slope index about —3. It can be found that both the number 
fraction and highest energy of accelerated particles decrease as {Obu) decreases. We 
have also tried different values of initial energies (not shown), and find that the 
acceleration efficiency d ecreases f or ele ctrons with higher initial energies, which is 
similar to the results of ( IBurgessl . 120061 ). 



143 



a, b c d ef g hijk 




236 233 240 242 244 246 243 250 SB.O B8.5 E9.0 

K /(c/w pi ) fi s l 

Figure 4.3: A typical electron trajectory analysis that shows acceleration by mul- 
tiple mirroring between ripples. The top left panel displays the trajectory of the 
representative electron in x-z plane, overlapped with contour of B z magnetic field 
where the gray-scale is the same as that in Figure 1; The top right panel shows the 
position of the electron in x coordinate as a function of time; The bottom left panel 
illustrates the energy of the representative electron E/E as a function of x; The 
bottom right panel shows the dependence of electron energy E/E on time. 
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Figure 4.4: The energy flux spectrum of electrons at VL c {t = 120 for different averaged 
shock normal angles. The red solid line is for the case with shock angle (9bu) — 90°, 
the blue dot dashed line and the black dashed line are for the cases with {0B n ) — 60° 
and 75°, respectively. 
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The effect of varying level of magnetic turbulence variance is examined in Figure 
14.51 We compare three cases with different turbulence variances a 2 = 0.1, 0.3, and 
0.5, respectively. At the end of simulations, the final energy spectra are similar 
at low energies, with significant variations in the spectra only at energies higher 
than 2 keV. It is found that the energy spectrum is hardened at high energies when 
the turbulence variance is largest, which indicates that the large-scale turbulence 
is more important for accelerating electrons to high energies. We argue that colli- 
sionless shocks that move through magnetic turbulence with significant power leads 
to efficient electron acceleration to high energies since the motion normal to the 
shock front is enhanced. The reason is that the meandering of field lines is en- 
hanced, which allows the electrons to have a better chance to travel though the 
shock multiple times. 

We examine the effect of different turbulence correlation lengths in Figure 14.61 
We compare three cases with different sizes of simulation box L z = 2000, 1024, and 
AOOc/ujpi, respectively. In each case, the correlation length is made to be the same 
as the size of the simulation box in ^-direction L c = L z . It is shown that for the 
case L c = L z = 2000c/a;pj, more electrons are accelerated to high energy. In the 
case of L c = L z = 400c/u pi , less electrons are accelerated and the highest energy is 
less than the other two cases. The more efficient acceleration for the case that L c 

e normal to the shock front 



is larger can also be understood as the motion of partic . 
is more enhanced. As derived by Jokipii and Parker (I1969I ). large-scale magnetic 
fluctuation dominate the random walk of magnetic field lines and can be expressed 
as 



(Ax 2 ) Pjk = 0) 
Az ~ Bl ' [4 - Z) 

where Ax and Az represent the displacement of a field line in the direction transverse 
and along the mean magnetic field, respectively. P(k = 0) represents spectral power 
at wave number k = 0. Although we use a restricted simulation box, our simulation 
is qualitatively consistent with this picture. In figure I4~7T1 we plot three magnetic field 
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Figure 4.5: The energy flux spectra of electrons at Qdt = 120 for an averaged 
perpendicular shock with different turbulence variances. The black dashed line, red 
solid line, and blue dot dashed line are in the cases that cr 2 = 0.1, 0.3, and 0.5, 
respectively. 



147 




Figure 4.6: The energy flux spectrum of electrons at Vt p it = 120 for averaged per- 
pendicular shock with different turbulence coherence lengths. The black dashed 
line, red solid line, and blue dot dashed line are in the cases that L c = 2000c/ 'cj pi , 
lOOOc/upi, and 400c/u pi , respectively. 
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lines generated from magnetic turbulence in the cases of three different coherence 
lengths. In this plot the horizontal axis x represents the spatial distance transverse 
to mean magnetic field with a plot range from — 1000c/ u P i to 1000c/ u P i, and the 
vertical axis z represents the spatial distance along the mean magnetic field. In 
order to show the difference between these field lines, the spatial range for z axis 
is changing for different field lines. In the three cases of L c = 14000c/ 'bj pi (blue 
dashed line), L c = 2000c/u;pj (red dotted line) and L c = 400c/ u pi (black solid line), 
the spatial ranges in z direction are made to be — 14000c/ u P i, — 2000c/u;pj and 
— 400c/uj p i, respectively. One can see that the x-direction displacement of the field 
line in the case of L c = 14000c/ u P i is much more enhanced compare to the case of 
L c = 400c/co> P j. Because the limited computational resource, the largest simulation 
is made to be L c = L z = 2000c/u P i, whereas for typical interplanetary parameters, 
this value is about lAOOOc/cjpi (0.01 AU). 

The idea of field line random walk (and its related large-scale irregular shock 
surface) is also useful for interpreting the observation of energetic particles associated 



with shock waves. An example is shown in Figure 14.81 which shows the profiles 
of the number of energetic electrons at Q C it = 100.0 as a function of x, for the 
case of (Obii) — 90°. The black solid line is the profile at z — 200c/ u P i, and the 
red dashed line shows the profile at z — 800c/u P i. The corresponding position of 
the shock front at each of these values of z are represented using dot lines. At 
z = 200c/u P i, it is observed that the energetic electrons travel far upstream up 
to about 100c/ uj p i from the shock. However, the profile at z = 800c/u P i shows no 
significant upstream energetic electron flux. The upstream energetic e lectron profiles 
show irregular features similar to in-situ observations reported by ISimnett et al. 
( 120051 ) (Figure 10). The irregular features are controlled by the global topology of the 
large-scale turbulent magnetic field lines, along which the accelerated electrons could 
travel far upstream. Additionally, energetic electron profiles in x direction generally 
show "spike-like" structure close to the shock front, which is usually observed in 
interplanetary shocks and Earth's bow shock. This feature is relatively stable within 
the simulation time once the upstream electron structure is developed. 
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Figure 4.7: Three turbulent magnetic field lines for different correlation lengths. 
The blue dashed line, red dot dashed line, and black solid line are in the cases that 
L c = 14000c/ Upi, 2000c/ Upi, and 400c/a;pj, respectively. 
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Figure 4.8: The profiles of the number of energetic electrons across the shock at 
z = 200c/u P i and z = 800c/ 'cu P i and time O ci t = 100.0, respectively. The red dot 
line and the black dot line label the correspoding positions of the shock fronts. 
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Since the temporal and spatial scales in our hybrid simulation is not large enough 
compared with the realistic scales of the magnetic turbulence in interplanetary space, 
we have to estimate how long the upstream electrons can be observed before the 
shock encounter in realistic parameters. We have carried out an order of magnitude 
calculation by assuming a planar perpendicular shock whose surface connects to 
meandering field lines of force at v arious places. The magn etic turbulence is assumed 
to be the two-component model (IMatthaeus et all Il990l ). discussed in Section I2~2l 
We take the value of coherence length in solar wind turbulence L c = 0.01 AU and 
the total variance of turbulence is 5B 2 /Bq = 0.3. One can estimate the diffusion 
of magnetic field lines transverse to the average direction of magnetic field D± = 
Ax 2 / Az. We take the value from numerical simulation that D i = 0.14L„ f Giacalone 



& Jokipii, Il999l ). Then we can derive that the field line wandering along the shock 
normal is Ax ~ 0.4L C . Using a shock speed of 800 km/s, we estimate that the 
electron foreshock region for a perpendicular shock propagating into a turbulent 
upstream is about 10 min utes. This roug h orde r of magnitude estimate agrees with 



Simnett et al 



(120051 ). The upstream ener getic electrons 



the in-situ observation by 
have been a lso observed by Voyager spacecraft, in a much lager scale ( iDecker et al. 



2005 



20081 ). Voyager II also observed the change of energetic electron intensity, 



presumbly due to large-scale magnetic field change. 



4.5.3 Comparison Between the Acceleration of Electrons and That of Ions and its 
Implication to SEP Events 

The correlation between electrons and ions in SEP events has been reported and 
discussed by a number of authors (see Section T4.2p . The correlation indicates that 
electrons and ions are accelerated close to the Sun by similar processes. We have 
shown that large-scale turbulence has important effects on accelerating electrons to 
high energy The result, along with the previous work on the acceleration of ions 
( IGiacalond . l2005a! Jbl). suggests that perpendicular shocks may play an important 



role in the acceleration of both electron and ions at shocks, which is consistent with 
the correlation between ions and electrons in solar energetic particle events (e.g., 
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Figure 4.9: The energy flux spectra of protons at Q C it = 120 for different averaged 
shock normal angle. The red solid line is in the case that the shock angle (#B n ) = 90°, 
the blue dot dashed line and the black dashed line are in the cases that (0Bn) — 60° 
and 75°, respectively. 
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Figure 4.10: The energy flux spectra of protons at Q C it = 120 for an averaged 
perpendicular shock with different turbulence variances. The black dashed line, red 
solid line, and blue dot dashed line are in the cases that cr 2 = 0.1, 0.3, and 0.5, 
respectively. 
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Figure 4.11: The energy flux spectrum of protons at Q p it = 120 for averaged per- 
pendicular shock with different turbulence correlation lengths. The black dashed 
line, red solid line, and blue dot dashed line are in the cases that L c = 2000c/ 'cj pi , 
lOOOc/ujpi, and 400c/u pi , respectively. 
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Oliver! . l20Q9h . 

Here we explicitly compare the energy spectra of electrons with that of protons. 
In Figure I4.9[ I4.10[ and I4.11[ we show the energy spectra of downstream protons 



from the hybrid simulations, corresponding to the spectra of electrons in Figure |4~4| 
14.51 and 14.61 We find that more efficient acceleration for protons can be obtained in 
the case of larger values of turbulence variances and correlation lengths. These agree 
well with the characteristics of the acceleration of electrons. In our simulations, the 
accelerated protons at the oblique shock with (9 Br,) = 60° are found to reach higher 



energi es, which is different from previous works (IJokipii 



1982 



1987 



Giacalone 



2005a 



bf). This is probably due to th e limite d temporal and spatial scales of our 



simulations. As shown by iGiacalond (l2005bl ) using test-particle simulations, the 
energy spectra of protons reach the highest energy in perpendicular shock case in 
a longer time scale Vt cl t ~ 50000 (this corresponds to 5 — 10 minutes for typical 
parameters in solar corona). However the current results from hybrid simulations 
do show a population of thermal protons can be accelerated to high energies in 
perpendicular shocks, which supports the idea that both electrons and protons can 
be efficiently accelerated by shocks with large shock normal angles. 



4.5.4 Summary 

We have presented the results of the acceleration of electrons (and also protons) 
at a perpendicular shock that propagates through a turbulent magnetic field. The 
acceleration of electrons are enhanced due to the effect of large-scale turbulence. 
The accompanying results for protons qualitatively show the correlation between 
accelerated electrons and accelerated ions in oblique shocks with large shock nor- 
mal angles. This indicates that quasi-perpendicular/perpendicular shocks play an 
important role in SEP events. 



This study will help to explain the correlation 



between electrons and ions in solar 



ener getic particles from both CME-driven shocks flCliverl. 120091) and solar flares f Shih 



et al.. 



20091 ). The result also poses a question on the contribution of parallel shocks 



in SEP events. Since it is difficult for parallel shocks to accelerate electrons, most 
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of SEPs are probably originated from perpendicular shocks. The a cceleration o: 



electr ons has to be included in current scenario of SEP events (e.g., iTylka et al. 
2005( 1 . 



4.6 The Effect of Large-Scale Magnetic Turbulence on the Acceleration of Elec- 
trons: Flare Termination Shocks 

In this section we use the combination of hybrid simulations and test-particle elec- 
tron simulation (see Section 14.41) to study the acceleration of electrons at flare ter- 
mination shocks predicted to exist in the vicinity of solar flares. Most of th is section 
has been published in the Astrophysical Journal ( Guo fc Giacalone . 20121 ). 



The existence of standing ter mination shocks has 
and numerical simulations (e.g. 



Shibata et al. 



1995 



jeen examine d by flare models 



Forbes 



1986 ). Solar fl a res ar e 



observed to be strong sources of energetic charged particles (lAschwandenl . 120021 ). 
The rele ase of magnetic energy by magnetic reconnection is thought to be the driving 
process (IMasuda et all Il994l ). While several mechanisms have been proposed to 



explain the acceleration o f charged particles in flares (see review by 



19971 : 



Zharkova et al 



Miller et al. 



2011 



and references therein), there is still no general consensus 
and this remains an unsolved problem. Recent hard X-ray observations of the non- 
thermal electron bremsstrahlung emission by Reu ven Ramaty High Energy Solar 
Spectroscopic Imager (RHESSI; 



Lin et al. 



20021 ) have provided more details of 



electron acceleration in solar flares. The observations indicate that a large fraction 
of released energy resides in high-energy electrons during a short amount of time. 
Hard X-rav sources above the top of magnetic loops have been detected fe.g.. Masuda 



et al., 



1991 



Krucker et al 



20101 ). providing important clue s to the acceleration 



Krucker et al 



(j2010h 



process. For example, the loop-top source recently reported by 
shows that a large number of electrons (> 5 x 10 35 ) are accelerated to more than 16 
keV and the highest energy reaches ~ MeV. Since the observed hard X-ray source 
requires very efficient acceleration, explaining how such a large number of electrons 
(probably also ions) are accelerated to high energy poses a challenge to theoretical 
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astrophysics. 

The existence of fast shocks in t 
dieted in flare models ( 



1986 



1988 



Shibata et al 



Shiota et al 



2003 



le re connect ion outflow region has b een pre- 



1995) and nume rical simulations (e.g.. 



Workman et aJj 



Forbes 



20111 ). Using MHD numerical sim- 



ulations, iForbesI (119861 ) studied the formation of a standing termination shock when 
a high-speed jet driven by reconnection encounters a closed magnetic loop. The 
geometry of the flare termination shocks can be represented by Figure 14.121 The 
high-speed jet created in the reconnection out-flow region collides with the top of 
the magnetic loop and produces a fast-mode, standing termination shock. The re- 
sulting flare termination shock has a unit normal to its surface that points nearly 
perpendicular to the magnetic field. This is a perpendicular shock (i.e., th e angle 
betwe en the upstream magnetic field and shock normal vector Qbu = 90°). [Forbes 



(119861 ) predicts the existence of this shock with a compres sion ratio of 2.0 a nd an 

(120 lib 



Workman et al. 



upstream Mach number as high as 2.3. A recent study by 
shows similar results. The observational evidence of the existence of flare shocks 
has been presented by 



Aurass et al 



(120021 ). 



Particle acceleration at flare termination shocks has been considered by a number 
of authors. It is usually thought that the injection of electrons is a problem at 
perpendicular shocks. We addressed this in previous sections. There has been some 
work s on facilitating particle acceleration at flare termination shocks. Tsuneta & 
Naito (119981) considered e l ectro n heating by a slow-shock pair as a pre-energization 



process. 



Somov fc Kosugil (119971 ) considered the role of plasma heating and collapsing 



magnetic trap at reconnecting magnetic field lines. In this paper we present results 
from a combination of hybrid simulations and test-particle electron simulations to 
study the electron energization at a flare termination shock in the existence of 
upstream magnetic fluctuations. We consider the nonlinear modification of a planar 
shock front by upstream Alfvenic fluctuations and its effect on electron acceleration. 
Although the plasma waves and turbulence in the reconnection outflow region could 
be considerably different from our simplified model, some intrinsic characteristics 
of this interaction, such as the braiding of magnetic field lines and shock rippling, 
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Figure 4.12: A cartoon illustration of the geometry of flare termination shocks. 

should still be preserved. We show that after considering a fluctuating upstream 
region, the electron acceleration in flare termination shock is a rapid and efficient 
process. A large fraction of the initial thermal electrons is accelerated to hundreds 
of keV and even reaches MeV energies in a very short time. This indicates that 
collisionless shocks may play an important role in particle acceleration in solar flares. 
In Section 14.6.11 we describe the initial conditions and parameters used in this paper. 
Section 14.6.21 discusses the simulation results. 

4.6.1 Initial Conditions and Parameters 

The size of the simulation domain L x x L z for each case is listed in Table 14.11 The 
flare termination shock is modeled by injecting plasma continuously from one end 
(x = 0) of the simulation box and colliding the reflecting boundary at the other end 
(x = L x ). The total variance of magnetic turbulence in each case is listed in Table 
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\. In all the cases, we take L c = L z . We expect large-scale magnetic turbulence 
to exist in reconnection outflow plasma. These fluctu ations can be t rigger ed by 



Krucker et al 



2008 



, for a 



reconnection, foot-point motion or other processes (see 
detailed discussion). Although this simplified form of magnetic-field fluctuations 
may not be realistic, any turbulence with large variances should allow strong field- 
line wandering, which is essential in our particle acceleration model. 

In order to produce a low Mach number shock as predicted by other numerical 
simulations, the inflow Alfven Mach number is taken to be M^o = 1-0. After 
reflection at the right boundary, this produces a shock with averaged Mach number 



of about 2.0 in the shock frame, consistent wit 



i the flare termination shock predicted 



by previous MHD simulations ( IForbesl . Il986l ; I Workman et al.l . 1201 ll ). The grid sizes 
are Ax = Az = 0.5c/uj p i and the time step is taken to be At = O.Olfi" 1 , where c/u P i 
is the ion inertial length and is the ion gyroperiod. The plasma beta and /? e 
are taken to be 0.03 and the ratio between light speed and upstream Alfven speed 
c/Vaq = 410, which roughly corresponds to an initial situation with temperature 
2 x 10 6 K, number density 8 x 10 9 cm" 3 and magn etic field Bq = 30 G, similar to 



constraints from observations (IKrucker et al. 



2010l ). Under these parameters, the 



average shock speed in the shock frame (also the jet outflow speed) is about 1460 
km/s. This speed is measurable because the outflow plasma from reconnection is 
moving at the Alfven speed. The estimate from observations (ITsuneta et all 119971 ) 
is roughly consistent with this value. 

In the second step, we integrate the relativistic equations of motion for an en- 
semble of test-particle electrons in the two-dimensional time-dependent electric and 
magnetic fields obtained in the hybrid simulations. These test-particle electrons are 
treated as a different part from the electron fluid in the hybrid simulations. We 
use a second-order spatial interpolation and linear temporal interpolation to get the 
field at the particle position. Initially we release a Maxwellian distribution with 
T e = 2.0 x 10 6 K in the upstream frame. The test-particle electrons are released 
upstream at Vt ci t = 30 after the shock has fully formed and far from the boundaries. 
The simulation domain in the x-direction is large enough so that no test-particle 
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Run 


L x {c/Upi) x L z {c/LOpi) 


dB /B 


"TO/ It? v i c U„~\A 
I/O \Fj > lo ke V J 


1 


500 x 400 


0.0 


1.3 


2 


500 x 400 


0.03 


4.5 


3 


500 x 400 


0.1 


8.0 


4 


500 x 400 


0.3 


9.8 


5 


500 x 800 


0.03 


4.9 


6 


500 x 800 


0.1 


8.9 


7 


500 x 800 


0.3 


11.9 



Table 4.1: Some parameters for different simulation runs. The size of the simulation 
domain, the variance of injected magnetic fluctuation, and the fraction of electrons 
whose energy is more than 15 keV at the end of simulation. 

electrons escape from the system. Strictly speaking, this test-particle simulation is 
only valid when the influence of the accelerated electron to the background fluid is 
negligible. However, in the end of the simulation, the initial Maxwellian distribu- 
tion of test-particle elections has been considerably changed due to the energization 
process at the shock front. This indicates that our approach may not be suitable in 
studying the long-term evolution of the termination shock. 

4.6.2 Simulation Results 

In this work we discuss the electron acceleration in flare termination shocks, which 
have low Mach numbers and high shock speeds. The plasma in solar corona is in a 
high temperature and low plasma beta (strong magnetic field) regime, which is con- 
siderably different from that in interplanetary space. We focus on the modification 
of shock surface by upstream Alfvenic fluctuations and its effect on acceleration of 
electrons. Table I47T1 lists some key parameters for all the simulation runs including 
the size of the simulation domain, the variance of the injected magnetic fluctua- 
tion, and the fraction of electrons whose energy is more than 15 keV at the end of 
simulation. For runs 1-4 we consider the effect of different variances of magnetic 
turbulence. The turbulence variances are from 0.0 to 0.3 and the sizes of the sim- 
ulation domain L x x L z = 500c/u P i x 400c/u P i (1.27 km x 1.02 km) for these four 
cases. For runs 5-7, the magnetic variances are the same as runs 2-4, but the size 
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of the simulation box is changed to L x x L z = 500c/u P i x 800c/u P i (1.27 km x 2.03 
km) to examine the effect of changing the coherence length. In the flare region, 
the strong large-scale Alfvenic magnetic fluctuation can be triggered by reconnec- 



tion, and cascade to small scales. This process is usually assumed to 
of magnetic turbulence requ ired in many acceleration models (e.g., 



3e the source 



Miller et al. 



1996 



Petrosian fc Liu. 



20041 ) . Since the size of our largest simulation domain is still 



much smaller than the observed hard X-ray emission region [L ~ 10 3 km), we do 
not consider the realistic geometry of flare termination shock but approximate it 
locally as a perpendicular shock that propagates into a plasma containing magnetic 
fluctuations. 

Figure HTT3] shows the color-coded contours of (a) magnetic field in the z direction 
B z and (b) ion number density n p from run 3 at VL ci t = 110.0. The magnetic field 
and plasma density have been normalized using the average upstream magnetic 
field £>io and in- flow density Uq. The averaged Alfven Mach number in the shock 
frame is about 2.0 and the average compression ratio i s abou t 2.1. As noted in 



the earlier works (IGiacalone . 



2005b 



Guo fc Giacalond . I201CH ) . the shock surface 



becomes distorted due to the interaction between the shock front and the upstream 
turbulence. Meandering magnetic field lines cross the shock front at various locations 
along the shock, which allows the electrons to cross and/or get reflected at the shock 
front multiple times. The shock-front rippling has also been s hown to contribute 
to pa rticle acceleration by mirroring electrons between ripples ( iGuo fc Giacalone . 



201fJ). 



Figure 14.141 presents the energy spectra dJ / dE of electrons. The green solid 
line shows the initial distribution of thermal electrons in the upstream region. The 
black solid line displays the energy distribution for all the electrons in downstream 
region at the end of simulation for run 1. In this case no pre-existing fluctuation 
is considered and the electron energization is primarily due to heating at shock 
layer. We have calculated that the effective electron kinetic temperature jump in 
the downstream region including the superthermal distribution is about 6 times of 
the upstream temperature. The simulated electron temperature jump is about 40% 



162 



400 



300 



200 



1 00 




350 400 450 

X /(c/W) 



400 



300 



200 



1 00 




350 400 450 

x /( c / w Pi) 



Figure 4.13: The color-coded images of (a) magnetic field B z , (b) ion density n p for 
run 3 at Q r d = 110.0. 



of the proton temperature jump across the shock layer in our hybrid simulation. 



This is consistent with the theoretica 



fast shocks is less than that of ions ([Goodrich Sz Scudder . 



prediction that the heating of elect r ons in 



19841 ; 



Scudder 



1223) 



and the observational constraints from measurements at planetary bow shocks and 



interplanetary shocks (jThomsen et all 119871 ; iSchwartz et all 119881 ). We note that 
in our simulation electron heating may not be determined very accurately since the 
test-particle electrons have no feedback to the electric and magnetic field at the 
shock layer. In the following we focus on the nonthermal acceleration of electrons at 
shocks after considering the pre-existing magnetic fluctuations. The blue solid, dot 
and dashed lines in Figure 14. 141 represent energy distribution for all the electrons in 
downstream region at the end of simulation (VL ci t = 130.0) for runs 2-4, respectively. 
At this time the energy spectra do not evolve anymore. It can be seen that the 
electrons are accelerated to high energy after considering the upstream magnetic 
turbulence. For higher variance of magnetic turbulence, there are more particles ac- 
celerated to high energy. For run 4 (5B 2 /Bq = 0.3), 9.8% of electrons are accelerated 
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Figure 4.14: The energy spectra of electrons at the end of the simulation (Qit = 
130.0). The energy spectra are normalized using N e v t h e {u P i/ c) 2 / keV , where N e is 
the total number of electrons used in the simulations and v t h e is the initial electron 
thermal speed. The green solid line shows the initial distribution of thermal electrons 
in the upstream region. The black solid line displays the energy distribution for all 
the electrons in downstream region at the end of simulation for run 1. The blue 
solid, dot and dashed lines represent results from runs 2, 3, and 4, respectively. 

to more than 15 keV at the end of the simulation. The efficient electron acceler- 
ation can be understood as stronger magnetic turbulence allows stronger field-line 
meandering, and the electrons can be taken by field lines of force and therefore gain 
energy at a hock front multiple times. 

In Figure 14.151 we examine the effect of changing the coherence length of the 
magnetic turbulence and focus on the high energy part of the energy spectra. It 
shows results from runs 5-7 (red lines, L z = 800c/ u p i) along with corresponding runs 
2-4 (blue lines, L z = 400c/ u pi ). It is shown that for larger coherence length, the 
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electrons could reach higher energy and the spectral slope tends to be flatter. The 
more efficient acceleration in runs 5-7 can be understood as the larger simulation 
domain in the direction of magnetic field allows more field line wandering normal 
to the shock (Ax 2 oc Az, where Ax is the field-line random walk normal to the 
averaged magnetic field and Az is distance along the field) therefore the electrons 
move across the shock more easily. This dependence shows that long-wavelength 
fluctuations are important to accelerate electrons to high energy. 

We also analyze the acceleration of protons, which are treated self-consistently 
in this problem (i.e., they are included in the hybrid simulation). Figure 14.171 shows 
the differential energy spectra of protons in the shock frame for runs 2-7. Similar to 
Figure 14.151 the results from runs 5-7 are represented by red lines and the results 
from runs 2-4 are displayed using blue lines. The accelerated protons show a similar 
dependence on turbulence variance and coherence length to that of electrons. This 
dependence has been found pr eviously fo r the ca se of higher Mach number shock 
and larger correlation length (jGiacalonej . l2005bl ). These results show that both 
electrons and protons can get efficiently accelerated. However, for the parameters 
we use the slopes of the energy spectra of protons are considerably steeper than 
that of the spectra of electrons. This is prob ably due to the lim ited temporal and 
spatial scales of our simulations. As shown by iGiacalond (j2005bl ) using test-particle 
simulations, the energy spectra of protons reach the highest energy in perpendicular 
shock case in a longer time scale tt&t ~ 50000. However the current results from 
hybrid simulations do show a population of thermal protons can be accelerated to 
high energies in perpendicular shocks, which supports the idea that both electrons 
and protons can be efficiently accelerated by shocks with large shock normal angles. 



4.7 Discussion and Conclusions 

In this chapter we studied the acceleration of electrons at collisionless shocks by 
utilizing a combination of a 2-D hybrid simulation to obtain the shock structure 
and a test-particle simulation to determine the motion of electrons. The hybrid 
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Figure 4.15: The energy spectra of electrons at the end of the simulation (fi^i = 
130.0). The energy spectra are normalized using N e v t he(^pi/c) 2 /keV, where N e is 
the total number of electrons used in the simulations and Vthe is the initial electron 
thermal speed. The red solid, dot and dashed lines represent the energy distributions 
for all the electrons in downstream region at the end of simulation for runs 5,6, and 
7 respectively. The blue solid, dot and dashed lines represent results from runs 2, 
3, and 4, respectively. 
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Figure 4.16: The relation between the turbulence amplitude 5B 2 /Bq injected in 
hybrid simulation and the percentage of electrons eventually accelerated to more 
than 15 keV. 
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Figure 4.17: The energy spectra of protons downstream of the shock at the end of 
the simulation {fl{t = 130.0), normalized using N p v t hp{upi/c) 2 /keV, where N p is the 
total number of protons used to plot the spectra and v t h P is the initial proton thermal 
speed. The red solid, dot and dashed lines represent the energy spectra for protons 
in downstream region at the end of simulation for runs 5, 6, and 7 respectively. The 
blue solid, dot and dashed lines represent results from runs 2, 3, and 4, respectively. 
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simulation provides realistic electric and magnetic fields within the transition layer 
of the shock that effect the motion of test-electrons, which is determined by solving 
the equation of motion. The interaction of the shock with pre-existing upstream 
fluctuations, and other nonlinear processes occuring in the hybrid simulation lead to 
a "rippling" of shock surface that also effects the transport of the electrons. We find 
that the electrons are efficiently accelerated by a nearly perpendicular shock. The 
turbulent magnetic field leads to field-line meandering that allows the electrons to 
cross the shock front many times. The rippling of the shock front also contributes 
to the acceleration by mirroring electrons between the ripples. This acceleration 
process is more efficient at perpendicular shocks. As (#B n ) decreases from 90° , 
both the number fraction and highest achievable energy of accelerated particles 
decreases. Based on our calculations, we conclude that perpendicular shocks are 
the most important for the acceleration of electrons. The current study is helpful in 
understanding the injection problem for electron acceleration by collisionless shocks. 
It is also found that different values of variances and correlation lengths of the 
injected magnetic turbulence also strongly affect the attainable maximum energy 
and accelerated fraction of electrons. The cases with larger turbulence variances 
or larger turbulence correlation lengths have a flatter energy spectrum than that 
with smaller turbulence variance, which suggests the enhanced motion of electrons 
normal to the shock front, due to enhanced field-line random walk, is of importance 
for the acceleration of electrons to high energies. 

Here we discuss the applications of this process to interplanetary shocks and 
SEP events (Section I4.7.ip and particle acceleration in solar flares (Section 14.7.2)) . 

4.7.1 Application to Interplanetary Shocks and SEP events 

We presented the results of electron acceleration at shocks under parameters similar 
to interplanetary shocks. The results can also be useful for considering acceleration 
of electrons by shocks in SEP events. For the case that the averaged shock normal 
angle (#e n ) = 90° and turbulence variance a 2 = 0.3, the electrons can be readily 
accelerated to up to 200 — 300 times the initial energy. The resulting spectrum is 
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flat between about 0.1 keV to 0.7 keV. At higher energies, the spectrum falls off 
with energy following a power law with a spectral slope of about —3. We found 
that the energetic electron density upstream and downstream of collisionless shocks 
show filamentary structures (Figure H~2~]) . This could help explain electron spike-like 
events observed upstream and downstream of terrestrial and interpla netary shocks 



( I Anderson et al. 



1979 



Tsurutani fc Lin 



1985 



Simnett et al. 



20051 ). Observation 



by Voyager 1 at the termination shock and in the heliosh eath also show the ey idence 
of electron spike-like enhancements at the shock front (IDecker et all 120051 ). The 
upstream spatial distribution of energetic electrons shows irregular features that 
depend on both the irregularity in the shock surface and the global topology of 
magnetic field lines. At first the electrons are accelerated and reflected at the shock 
front, and then they travel upstream along the magnetic field lines. The electrons 
could be taken far upstream by field line random walk. This result can possibly lead 
to an interpretation to the complex electr on foreshock events recently obse r ved t o 



be associated wit h int erplanetary shocks (|Bale et al. 



Bale et al 



1999 



Pulupa fc Bale . 



2008 . 



( 119991 ) and iPulupa fc Bald (120081 ) proposed that the complex upstream 
electron events result from large-scale irregularities in shock surface. In this paper 
we have demonstrated that the upstream electron flux may be controlled by both 
an irregular shock surface and by large-scale meandering magnetic field lines. 

We have also presented the accompanying results for protons that qualitatively 
show the correlation between accelerated electrons and accelerated ions in oblique 
shocks with large shock normal angles. The shocks can efficiently accelerate both 
electrons and ions. This indicates that quasi-perpendicular/perpendicular shocks 
play an important role in SEP events. 



4.7.2 Application to Particle Acceleration in Solar Flares 

Understanding particle acceleration in solar flares is a challenge since only remote 
observations are available and it is hard to identify the main mechanism. While 
it is commonly thought that magnetic reconnection drives the energy release, the 
detailed physical process involved in accelerating the electrons and ions is still not 
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clear. We studied electron acceleration at flare termination shocks that have been 
predicted by numerical simulations and flare models. We find that electrons are 
rapidly and efficiently accelerated at a flare termination shock in the presence of pre- 
existing magnetic fluctuations. Electrons are accelerated to a few MeV in 100 ion 
gyroperiod (of the order of a millisecond) and more than 10% of thermal electrons are 
accelerated to more than 15 keV given a sufficiently strong magnetic turbulence. We 
also show that electron acceleration is more efficient for larger turbulence variance 
5B 2 1 Bq and/or a larger turbulence coherence length L c . Both of these indicate that 
large-scale field-line meandering plays an essential role in accelerating electrons at 
a shock front. Our simulations show that after considering the magnetic turbulence 
the flare termination shock could ac celerate elect r ons t o much higher energies than 
usual drift shock acceleration (e.g., iMann et al.l . 120091 ). We note that the similar 



mechanism has been shown to efficiently accelerate ions and has similar dependence 
on the turbulence properties. This correlation between ions and electrons is actually 
commonly observed in solar energetic particle events. For the parameters used 
in our simulations, the accelerated protons have energy spectra steeper than that 



of electrons. This is different 



interplanetary space (IGiacalone 



rom the previous results for p arameters similar to 



2005b 



Guo fc Giacalonei . 120101 ) . We note that these 



results are carried out for energies lower than the injection energy for diffusive shock 
acceleration. When the pitch-angles of charged particles are scattered sufficiently 
as to be trapped near the shock, the energy spectra of the accelerated particles are 
presumably close to that predicted by diffusive shock acceleration. 

We also note that for the situation we study, the resulting distribution of elec- 
trons is non-Maxwellian. The structure of collisionless shocks may be considerably 
modified by accelerated particles. While this effect is not considered in our test- 
particle simulations, it may be important. Due to the limitation of computation, a 
full particle simulation of collisionless shock with realistic mass ratio mj/m e = 1836 
that includes the influence of turbulent upstream magnetic field is not available so 
far. The evolution of this flare termination shock remains to be explored. Also, 
other plasma effects like emission process in flare region may need to be considered 
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to directly compare with the observations. 
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CHAPTER 5 
Conclusions and Future Work 



5.1 Summary and Conclusions for the Dissertation 

In this dissertation, we studied the effects of magnetic fluctuations on the accelera- 
tion and transport of charged particles in the heliosphere. We started by reviewing 
the basics of the acceleration and transport of charged particles and r elevant physic s 



in the heliosphere. We discussed the Parker's transport equation (jParkerl . Il965l ). 
which has been the fundamental equation to study the transport and acceleration 
of energetic charged particles. Then we focused on the limitation of the transport 
equation and discussed the possible solutions. Since Parker's transport equation as- 
sumes a quasi-isotropic distribution function in momentum space, in the case that 
the distribution function of energetic particles is highly anisotropic, the evolution of 
the distribution function cannot be properly described by the transport equation. 
The examples are the initial release of solar energetic particles and the acceleration 
of low-energy particles at the shock front. In addition, the observations of impulsive 
SEP events show fine structures in intensity-time plots on small temporal scales 
(hours) that cannot be easily described by a large-scale spatial diffusion. 

In Chapter 2 we presented numerical simulations for the propagation of SEPs 
in the inner heliosphere. We numerically integrated the trajectories of energetic 
charged particles in the turbulent magnetic field generated from two commonly 
used magnetic turbulence models (the foot-point random motion model and the 
two-component model). The observations of SEP events are simulated by collecting 
charged particles which reached 1 AU. We study the velocity dispersion of SEPs in 
the turbulent magnetic field and estimate the error involved in the onset analysis. 
We find that the the velocity dispersion can be well produced by this model. For 
a typical turbulence variation 5B 2 /Bq ~ 0.1 observed at 1 AU and a large source 
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region, we find that the difference between the apparent release time inferred from 
the onset analyses and the actual release time is less than a few minutes, but the 
apparent path length can be significant different than the real path length along 
the average magnetic field line. For the foot-point random motion model, the error 
for the inferred release time is smaller than that of the two-component model. We 
have also reproduced SEP dropouts in the numerical simulations using the foot 
point random motion model, assuming the SEP source region is smaller than the 
correlation scale. The widths of these dropout are typically several hours, similar 
to the time scales of dropouts observed in space. The velocity dispersion of the 
energetic particles appears to have different path lengths, which indicates that the 
energetic particles travel along different field lines. We have also attempted to 
use the two-component model to numerically simulate the dropouts of energetic 
particles. However, we rarely find the evidence of SEP dropouts in our simulation 
for the two-component model. This is because the parallel diffusion coefficient of 
particles in the two-component model is considerably smaller than that in the foot- 
point random model. This result questions the popular used two-component model 
in that it gives more pitch-angle scattering than that constrained by the observation 
of SEP dropouts. 

In Chapter 3 we studied two processes for particle acceleration at shock waves. 
The first problem is associated with the acceleration of low-energy particles at 
shocks. We presented a numerical study on the acceleration of thermal protons 
at parallel shocks using 3-D hybrid simulations. The 3-D simulations removed the 
artificial restriction of the motion of charged particles in previous 1-D and 2-D sim- 
ulations. The results confirmed the injection mechanism at parallel shocks that the 
accelerated particles are originated from reflected particles. In the second study we 
illustrate the effect of a large-scale sinusoidal magnetic field variation. This simple 
model allows a detailed examination of the physical effects. As the magnetic field 
lines pass through the shock, the connection points between field lines on the shock 
surface will move accordingly. We find that the region where connection points ap- 
proaching each other will trap and preferentially accelerate particles to high energies 
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and form "hot spots" along the shock surface. The shock acceleration will be sup- 
pressed at places where the connection points move apart each other. Some of the 
particles injected in those regions will transport to the "hot spots" and get further 
accelerated. The resulting distribution function is highly spatial dependent at the 
energies we studied, which could give a possible explanation to the Voyager obser- 
vation of anomalous cosmic rays (ACRs). This mechanism gives an interpretation 
for the observation that the ACRs did not saturate at the termination shock. 

In Chapter 4 we studied the acceleration of electrons at collisionless shocks by 
utilizing a combination of a 2-D hybrid simulation to obtain the shock structure 
and a test-particle simulation to determine the motion of electrons. We find that 
the electrons are efficiently accelerated by a nearly perpendicular shock when the 
large-scale pre-existing magnetic fluctuations are considered. The turbulent mag- 
netic field leads to field-line meandering that allows the electrons to cross the shock 
front many times. The rippling of the shock front also contributes to the accelera- 
tion by mirroring electrons between the ripples. This acceleration process is more 
efficient at perpendicular shocks. As (9bu) decreases from 90°, both the number 
fraction and highest achievable energy of accelerated particles decreases. Based on 
our calculations, we conclude that perpendicular shocks are the most important for 
the acceleration of electrons. The current study is helpful in understanding the 
injection problem for electron acceleration by collisionless shocks. It is also found 
that different values of variances and correlation lengths of the injected magnetic 
turbulence also strongly affect the attainable maximum energy and accelerated frac- 
tion of electrons. The cases with larger turbulence variances or larger turbulence 
correlation lengths have flatter energy spectra than the cases with smaller turbu- 
lence variances, which suggests that the enhanced motion of electrons normal to the 
shock front, due to enhanced field-line random walk, is of importance for the accel- 
eration of electrons to high energies. We discussed the applications of this process 
to interplanetary shocks and SEP events (Section 14. 7. ip and particle acceleration in 
solar flares (Section |4.7.2p . We also discussed the implication of this study to solar 
energetic particles (SEPs) by comparing the acceleration of electrons with that of 
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protons. The intensity correlation of electrons and ions in SEP events indicates that 
perpendicular or quasi-perpendicular shocks play an important role in accelerating 
charged particles. 



5.2 Future Work 



5.2.1 Effect of Shock Geometry on the Acceleration of Charged Particles in Gradual 
SEP Events 

The acceleration of solar energetic particles (SEPs) remains to be one of the most 
important unsolved problems in heliospheric physics. Observations indicate that the 
acceleration at high energies is highly variable in spectral properties and elemental 
composition. Many recent works suggest that effects like seed particles, shock ge- 
ometries and/or the generation of self-exci ted waves can play an important role on 
the observed variable energy spectra (e.g., iTylka et al.l . 120051 ). 

In Section 3.5 we have considered the effect of large-scale shock geometry on 
the diffusive shock acceleration (DSA) of anomalous cosmic rays at the termination 
shock. The results show that this effect can significantly modify the well-known 1-D 
steady state solution of DSA. Similarly, we plan to study numerically the acceleration 
of particles in CME-driven shocks, including the effect of shock geometry. Since the 
CME-driven shocks are known to have large-scale non-planar shapes, we expect that 
the particles can sample shock fronts with various shock normal angles during the 
acceleration. This effect has not been considered by previous works, which usually 
assume a planar shock or no cross-field d iffusion (e.g. 



Sandroos fc Vainio 



2007 



Li et al. 



Lee 



2005 



Tylka et al. 



2005 



20091 1 . This work will treat particle acceleration 



in CME-driven shocks in a more realistic way. 



5.2.2 Understanding the Physical Processes in the Acceleration of Low-Energy 
Particles at Shocks 

Solving the acceleration of energetic particles in the heliosphere requires the knowl- 
edge of the acceleration of low-energy particles. Since the DSA is not concerned with 



176 



low-energy particles that has high anisotropies, most of DSA models do not include 



the acceleration o 



1981 



Ellison et al. 



Ellison 



low-e nergy particles or treat it in a ad hoc way (e.g 
Il990 ). To fully understand the observed energetic particles, we 
need to model the acceleration of low-energy protons, electrons, and also heavy ions. 

a) Recent observations have indicated that the pre-accelerated particles from 
solar flares or other source may be important in acceleration of SEPs. In Chap- 
ter 3 we have used 3-D hybrid simulations to study the acceleration of low-energy 
protons at collisionless shocks. We plan to further explore the injection problem 
using the self-consistent hybrid simulations including the seed particle population 
pre-accelerated in solar flares or other sources. 

b) We also plan to use full particle simulations (kinetic ions and kinetic electrons) 
to study the acceleration of electrons at collisionless shocks. In Chapter 4 we have 
shown that low-energy electrons can be efficiently accelerated when a shock prop- 
agates normal to a magnetic field that contains pre-existing large-scale magnetic 
turbulence. The simulation model include a combination of a hybrid model (kinetic 
ions and fluid electron) and a test-particle electron model. Since the model does not 
consider the electron-scale plasma physics and the feedback of accelerated electrons 
on the shock front, we need to use a full particle simulation model to consider the 
effects in order to fully solve the problem of acceleration of electrons at collisionless 
shocks. This work will consider the interplay between large-scale effects (large-scale 
magnetic fluctuations and shock ripples) and small-scale effects (whistler waves at 
shock fronts). 

5.2.3 The Application of the Hybrid Simulations in the Heliospheric Plasma Pro- 
cesses 



During the doctoral study, we have improved the parallelization of the 



3-D hybrid simulation models f lGiacalone et al. 



2000 



iGiacalone 



2004 



-D, 2-D and 



2005b|). The 



new version of the code has good scalability and has been tested on NASA's Pleiades 
supercomputer using a few thousand CPU cores. The high performance of the code 
on supercomputers allows us to study a variety of problems using multi-dimensional 
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simulations or 1-D simulations with large simulation boxes and long simulation 
times. While we have presented some results using the new code, plenty of other 
plasma physics problems in the heliosphere can be studied using this numerical tool. 
We h ave plans to study num erically the ion cyclotron waves excited by fresh pickup 



ions (IFlorinski et al. 



20101 ) using 2-D hybrid simulations, plasma instabilities and 



waves upstream and downstream o f shocks, and the acceleration of charged particles 
by multiple shocks (ILi et al.l . |2012| ). 
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